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INTRODUCTION 

Studies  have  shown  thd  individuals  with  a  small  moment  of  inertia  (slender  bones)  are  at  higher  risk  of 
developing  stress  fractures  during  basictraining  (Beck  et  al.,  1996;  Milgrom  et  al.,  1989).  Becausethe 
size,  shape  and  material  properties  of  adult  bones  are  determined  early  in  life,  it  will  be  important  to 
understand  the  biological  processes  (or  factors)  that  contributeto  variation  in  these  morphometric  and 
material  traits  in  the  context  of  the  development  of  stress  fractures  Understanding  how  the  genetic 
background  affects  fracture  healing  will  I ead  to  further  ina'ghts  into  skeletal  tissue  response  to 
pharmacological  agents  that  might  promote  fracture  healing  and  predict  if  different  individuals  will  heed 
at  differing  rates.  The  hypothesis  of  these  studies  is  that  the  genomic  processes  that  lead  to  the  variations 
in  both  structural  and  material  properties  of  bone  development  will  be  recapitulated  in  the  developmental 
mechanism(s)  that  controls  the  bone’s  structural  geometry  and  materia  properties  during  fracture  healing. 
To  test  this  hypothesis,  si  mpl  e  transverse  fractures  were  generated  in  the  femora  of  three  strains  of  mice, 
A/J,  B6,  and  C3H.  These  stransof  mice  were  chosen  3  nee  they  have  known  differencesin  peak  bone 
mass  which  has  been  related  to  both  structural  and  materia  varidions  intrinsic  to  these  strai  ns  The 
ultimate  god  of  this  proposal  was  to  identify  the  genetic  components  that  regulatethe  variations  in 
development  of  bone  mass  and  geometry  through  the  comparison  of  fracture  repair  processes  in  the  three 
strains  of  mice.  Two  goals  were  proposed  to  address  the  research  hypothesis.  The  objective  of  goal  1 
was  to:  a)  Determine  the  nature  of  the  quantifiable  differences  in  the  biomechanical  and 
microarchitectural  properties  of  the  fracture  cal  I  uses  over  a  42  day  healing  period  in  the  three  strains  of 
mice  b)  Determine  how  the  genetic  strain  variations  in  the  microarchitectural  properties  relates  to  the 
rates  of  biomechanical  healing  as  well  asthe  relationship  between  micro  architecture  aid  regain  of 
biomechanical  strength.  The  objective  of  goal  2isto:  a)  Determine  the  underlying  genomic  processesthat 
are  activated  and  accompany  fracture  healing  in  the  three  genetic  strains  of  mice  and  identify  strain 
variation  in  the  patterns  of  mRNA  expression  by  carrying  out  both  select  and  full  transcriptional  profiling 
of  mRNA  expression  over  the  sane  healing  period  used  in  aims  one  b)  Identify  tissue  based  differences 
i  n  the  mechani  sms  and  assod  ated  genes  thd  ae  responsi  bl  e  f or  geneti  c  strai  n  vari  ati  ons  i  n  bone 
devei  opment  and  fracture  heal  i  ng. 

BODY 

Rrogresstowads  completion  of  Goal  1;  Mechanical  testing  and  micro  CT  has  been  completed  for  studies 
exami  ni  ng  al  I  three  strai  ns  of  mi  ce.  We  revi  sed  the  end  poi  nt  from  42  days  to  35  days,  si  nee  histologies 
assessments  of  the  42-day  post  fracture  calluses  showed  that  most  of  the  remodeling  in  the  mouse  calluses 
had  been  completed  and  it  was  decided  to  initially  obtain  mechanical  and  CT  properties  that  would  be 
more  reflective  of  the  period  of  primary  bone  remodeling  and  regain  in  mechanical  property  versus 
terminal  rates  of  regain  of  strength.  A  series  of  experiments  assessing  the  biomechanical  and  structural 
properties  of  healing  fracture  calluses  in  the  three  strains  of  mice  described  above  (A/J,  C57/B6,  and 
C3H)  were  carried  out.  Two  times  during  fracture  healing  were  initially  chosen  for  examination.  Day21 
was  chosen  since  it  is  one  of  the  earliest  times  at  which  mechanical  property  can  be  measured  and  is  a  key 
transition  point  between  the  endochondral  and  primary  bone  formations  peri ods  Day  35  was  examined 
since  it  is  representative  of  when  the  bone  should  have  regained  its  pre  fracture  strength  and  is  a  period  of 
intense  secondary  formation  and  remodeling.  MicroCT  was  carried  out  on  between  10-15  cal  I  uses  per 
drain  per  time  point  and  representdive  mi  croCT  images  of  the  three  strai  ns  at  the  two  time  poi  ntsae 
presented  (Figure  1).  MicroCT  images  show  that  at  21  days  post-fracture,  there  was  significant  vaiation 
i  n  the  size  and  structure  of  the  cal  I  uses  among  the  three  i  nbred  mouse  strai  ns.  The  average  total  a-ea  aid 
the  polar  moment  of  inertia  as  quantified  from  the  microCT  images  of  the  21  and  35  days  cal  I  uses  in 
comparison  to  the  intact  bones  for  A/J,  B6,  and  C3H  strains  are  seen  in  the  lower  panels  of  Rg  6. 

S  gnif  icant  differences  i  n  the  size  of  the  cal  I  uses  were  observed  among  the  i  nbred  mouse  strai  ns  d  21 
days  The  area  and  pda  moment  of  i  nerti  a  of  the  fracture  cd  I  us  for  the  A/J  mi  ce  wae  321%  and  1 022% 
greda  compaed  to  the  intact  bone.  I  n  contrast  the  ae a  and  poi  a  moment  of  inertia  of  the  B6  mouse 
fracture  cd  1  uses  wae  525%  and  2603%  greata  compaed  to  the  intact  bone.  The  aea  aid  pda  moment 
of  i  nerti  a  for  C3H  mice  fracture  cd  I  uses  wae  484%  and  2286%  greda  compaed  to  the  i  ntact  bone. 
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Even  more  striking  however  was  the  variation  (all  statistically  significant)  that  was  seen  when  these 
values  were  compared  between  mouse  strains.  The  microCT  images  also  showed  several  striking 
di ff erences  i  n  the  qual  ity  of  the  ti  ssue  that  w as  formed.  The  most  obvi ous  di ff erence  was  seen  i  n  the  I  ack 
of  bony  bridging  as  exemplified  by  the  area  of  X-ray  transiucency  in  the  central  region  of  the  C3H 
calluses  at  21  days,  which  would  be  indicative  of  cartilage  or  fibroustissue  still  being  present.  On  the 
other  hand  by  35  days  all  three  of  the  strains  of  mice  showed  bony  bridging  with  the  C3H  mice  showing 
thedensest  bone.  Torsional  teststofailurequantitatively  validatedthe  visual  observations  demonstrating 
that  the  C3H  femurs  had  a  stdistically  lower  mechanical  strength  and  stiffness  at  21  daysthan  either  the 
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Fig  l.  Structural  and  Biomechanical 
Compari  sons  of  Fracture  Cal  I  uses  From 
Three  I  n  Bred  Strains  of  Mice. 

Upper  Panels.  Representative  microCT 
images  in  both  transverse  and  longitudinal 
orientations  taken  at  21  and  35  days  post 
fracture.  The  orientation  of  the  microCT 
images  is  denoted  for  the  left  most 
transverse secti onoftheAJstrain.  Boxed 
area  of  C3H  day  21  i  mage  shows  area  that 
remained  unbridged  by  bone. 

Lower  Panels.  Biomechanical  and 


Structural  Data.  The  data  depicted  in  each 
of  the  four  paid  s  i  s  denoted. 

2  Wav  Anova  For  Genotype  and  Aoe 

Moment  of  Inertia  p=.0001  p=.0001 

Area  p=.0047  p=.0002 

Stiffness  p=.0001  p=. 0002 

Torque  to  Failure  p=.0194  p=.00Q2 
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Relative  Exf>re*sio 


A/J  or  B6  mice.  Interestingly,  the  B6  femurs  regained  their  pre-fracture  strength  by  21  da/s,  which  was 
due  to  both  to  the  greater  cross  sectional  areas  of  the  cal  I  us  and  the  tissue  transition  to  mineralized  bone. 
These  data  showed  then  that  the  mouse  with  the  strongest  pre  fractire  bones  (C3H)  initially  healed  the 
slowest  while  the  B6  mi  ce  heal  ed  the  fastest.  They  also  clearly  show  that  almost  all  of  the  structural 
features  of  these  three  strains  of  mice  are  genetically  inherited  and  recapitulated  during  the  processes  of 
fracture  healing.  More  importantly  they  show  that  there  will  be  genetic  variations  in  rates  of  bone 
healing.  This  will  have  eventual  clinical  importance  since  as  we  define  these  differences  in  animals  this 
will  lead  to  the  development  of  predictive  indices  to  identify  human  characteristics  that  will  be  diagnostic 
to  variable  healing  after  skeletal  trauma. 


Progress  towards  completion  of  Goal  2: 

The  objective  of  goal  2  were  to:  si)  Determine  the  underlying  genomic  processes  that  are  activated  and 
accompany  fracture  healing  in  the  three  genetic  strains  of  mice  and  identify  strain  variation  in  the  patterns 
of  mRNA  expression  by  carrying  out  both  select  and  full  transcriptional  profiling  of  mRNA  expression 
over  the  same  healing  period  used  in  a  ms  one  b)  Identify  tissue  based  differences  in  the  mechanisms  and 
associated  genes  that  are  responsible  for  genetic  strain  variations  in  bone  development  and  fracture 
healing.  The  experiments  that  refocused  on  tissue  based  processes  that  are  related  to  the  genetic 
variations  in  bone  healing  are  first  discussed. 
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Fig  2.  Cartilage  and  Bone  Extracellular  Matrix  mRNA  expression  Profiles  Across  the  Time  Course  of 
Fracture  Healing.  Panel  A)  Representative  RPA  analysis  of  osteoblast  and  chondrocyte  ECM  mRNA 
expression.  The  position  of  the  protected  bands  for  each  ECM  mRNA  is  denoted  in  the  figure  Panel  B) 
Overall  anabolic  activity  was  estimated  based  on  direct  beta  counting  of  all  of  the  expressed  ECM  bands 
normalized  to  the  L32  ribosomal  protein  housekeeping  gene,  while  the  percent  cartilagewasestmated  asthe 
ratio  of  thetotal  counts  of  CollOal  and  Col2a1.  Panel  C)  Temproal  profiles  of  expression  for  the  expression 
of  individual  mRNAs.  Each  mRNA  is  denoted  in  the  figure.  The  Bottom  two  panels  show  the  ratio  of 
Col10a1:Col2a1  aid  the  Col2a1:Osteocalcin  expression. 


6 


In  order  to  further  define  how  the  underlying  genetic  difference  of  the  three  mouse  strains  effects  the 
cellular  and  molecular  mechanisms  by  which  fractures  heal,  we  assayed  theecpression  of  a  group  of 
extracellular  matrix  mRNAsthat  would  be  indicative  of  either  cartilage  or  bone  differentiation  aid 
activity.  These  results ae  presented  in  Figure  Z  Both  a  qualitdive  and  quantity ve  assessment  of  the 
mRNA  expression  data  revealed  a  number  of  unique  features  of  theextracellula  marix  mRNA 
express  on  patterns  that  ae  i  nf  ormati  ve  t  o  the  cel  I  ul  a  and  biological  mechani  sms  that  ae  vari  abl  e  i  n  the 
three  strains  of  mice.  Thefirst  and  most  straight  forwad  observation  is  the  overall  diffaencesin  aiabolic 
activities  between  the  three  strains  of  mice  such  that  C3H>B6>AJ.  As  described  above  this  difference 
was  previously  shown  in  studies  of  Sheng  et  al.  (2004)  inthe  activity  of  osteoblasts  from  miureaiimds 
In  the  context  of  the  overall  types  of  anabolic  activity  it  is  also  intaesting  to  note  that  the  percentage  of 
cartilage  mRNA  expression  was  highest  in  the  AJ  strain  and  lowest  in  theC3H  strain.  This  would  suggest 
that  in  both  the  AJ  and  B6  drains  that  therewere  greeter  numbers  of  mesenchymal  stem  cellsthat  commit 
to  chondrocytes  and  progress  through  ai  endochondral  pathw^  of  bone  formkion  than  to  an  osteoblastic 
or  intramembraneous  pathway.  This  interpreted  on  isalso  consistent  when  the  ratio  of  the  terming 
orteoblart  marker  osteocalcin  was  calculated  relative  to  the  expression  of  the  cartilage  Col2a1  ma^er 
(lower  right  panel).  Theratioof  CollOal  toCol2aisalsoquiteinformativetothecellular  processesthat 
are  contri  buti  ng  to  f  racture  heal  i  ng  si  nee  it  i  s  i  nefi  cati  ve  of  the  rate  of  progressi  on  of  a  chondrocyte 
through  endochondral  differentiation.  In  this  regard,  one  of  the  most  interesting  aspects  of  these  data  was 
the  decreased  ratio  of  CollOal  to  Col2a1  in  theB6  call  uses  preceding  a  very  quick  rise  at  day  14.  TheA/J 
cal  I  uses  also  showed  a  peak  in  this  rati  oat  day  14,  suggesting  that  the  endochondral  process  is  completed 
soon  after  d^  14.  In  contra^,  theC3H  calluses  showed  a  very  gradual  and  continuing  increase  in  the 
CollOal  to  Col2a1  ratio  during  the  time  course  that  was  examined  in  thisstudy  demonstrating  that  in  this 
st rai  n  of  mi  ce  that  the  endochondral  process  takes  much  I  ortger  to  compl  eta 

Taking  ail  three  sets  of  data  together  shows  that  there  i  s  a  very  I  age  di  chotomy  of  geneti  c 
variation  in  the  biological  mechanisms  by  which  fracture  healing  proceeds.  I  n  the  case  of  the  most  rapid 
rate  of  healing  seen  in  B6  strain,  the  process  of  heeling  was  characterized  by  a  very  large  expaision  of 
immature  chondrocytes  as  exemplified  by  the  high  ratio  of  cartilage  to  bone  and  low  riio  CollOal  to 
C°l2a1,  followed  by  a  very  rapid  rate  of  hypertrophic  differentiation.  In  contrast  in  very  slowly  heating 
C3H  mi  ce  the  rati  o  of  carti  I  age  to  bone  was  the  I  owest  and  rate  of  hy pertrophi  c  maturati  on  was  the 
slowest.  In  summary  then  rapid  healing  was  achieved  by  having  a  robust  endochondral  mechanism  that 
woa  manife&ed  structurally  in  a  large  cross  sectional  callus  area  and  the  rapidity  of  progression  through 
the  endochondral  phase.  In  contrast  the  slowest  healing  C3H  mice  had  the  least  endochondral 
development  and  only  slcwly  progressed  through  the  endochondral  phase,  seen  grossly  in  the  delay  in 
bony  bndging  and  very  low  stiff  ness  of  the  tissue  at  21  days.  Even  though  the  cross  sectional  of  area  of 
the  C3H  mi  ce  was  ajbstanti ally  greater  than  that  of  the  A/J  stra  n  the  i  ncreased  cross  secti  onal  area  of  the 
cal  I  us  tissue  failed  to  achieve  increased  strength  and  thus  healed  slower  la-gely  because  of  the  poorer 
materi  ai  property  of  the  cal  I  us  due  to  the  per  si  stence  of  the  endochondral  carti  I  age  and  the  more  rapi  d 
progression  of  the  endochondral  process  in  the  A/J  and  B6  strains  Examination  of  14  day  post  birth 
growth  plates  demondrated  si  mi  la-  temporal  patterns  of  endochondral  progression  with  the  AJ  aid  B6 
mouse  strains  showing  ajbstantially  advanced  stages  of  maturation  as  compared  to  C3H  that  lajged 
behind.  These  data  suggested  thi  the  genetic  elements  which  control  the  primary  structural  attributes  of 
bone  quality,  are  recapitulated  during  fracture  healing  andean  be  easily  3SS6SS6d  in  the  temporal  window 
Of  fracture  healing.  These  data  indicate  that  the  basic  temporal  and  quantitative  attributes  of 
endochondral  bone  formation  during  fracture  healing  are  genetically  regulated.  Finally,  we  show  that  the 
rate  of  bone  healing  is  not  simply  related  to  variability  in  bone  quality  but  isalso  controlled  by  other 
factors  associated  with  geretic  variability  intrinsicto  endochondral  bone  formation. 
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These  data  drive  home  the  conclusion  that  the  much  slower  healing  of  the  C3H  mice  was  due  to  intrinsic 
genetic  variability  in  biological  mechanisms  that  effected  the  healing.  They  al  so  suggest  that  a  number  of 
the  underlying  mechanisms  controlling  the  genetic  variations  of  bone  quality  may  be  observable  in  the 
biological  processes  of  fracture  healing.  These  results  have  been  presented  in  three  abstracts  included  as 
appended  documents  1-3.  We  ha/e  also  used  funding  from  this  contract  to  partially  support  work  reported 
in  the  Appended  manuscript  6.  These  studies  establish  the  basic  methodology  to  collect  histologies* 
assessments  from  fracture  call  uses  as  well  as  demonstrate  thd  callus  formation  is  asymmetrical  to  the 
fracture  and  unique  to  developmental  growth  the  bone  that  was  fractured.  Thus  a  femur  which  grows 
from  the  di  Sal  epiphyi  si  swill  ha/e  a  cal  I  us  that  forms  asymmetrically  with  moS  of  the  endochondral 
bone  formation  on  the  distal  side.  In  contrast  thetibiathatgrowsfromitsprcximal  end  has  cal  I  us  that 
form  endochondral  bone  proximal  to  the  fracture. 

A  major  focus  of  goal  2  wasto  carrying  out  large  scale  transcriptional  profiling  comparisons 
between  the  three  strai ns  A  brief  summary  of  our  initial  studies  establishing  transcriptional  profiling 
analysis  across  fracture  healing  in  the  C57/B6  stra'n  are  presented  in  appendix  5.  The  micro  array 
experiment  of  the  comparison  of  the  three  strains  has  been  technically  completed,  but  data  analysis  is 
proceeding  now  and  will  continue  after  the  contract  isfirished  throughout  much  of  2006) 

Key  Research  Accomplishments 

1)  The  data  validate  our  hypothesis  that  the  genetic  variability  that  leads  to  the  differences  in  the 
quantifiabletraitsthat  effect  bone  quality,  also  lead  to  specific  quantifiable  differences  in  the  mechanisms 
by  which  a  fractured  bone  heals 

2)  The  data  identify  two  genetic  extremes  of  a  slow  and  fast  fracture  healer,  that  can  be  used  in  our 
studies 

3)  The  data  demonstrate  that  assay  of  complex  patterns  of  mRNA  expression  has  the  sensitivity  to  detect 
and  define  variations  in  biological  processes  that  are  genetically  inheri  able  (as  an  example,  variationsin 
patterns  of  endochondral  progression). 

Reportable  Outcomes 

Abstracts 

Jepsen  KJ,  Price  C,  Nasser  P,  A I  apart  M ,  Stapleton  SN,  Slkman  L,  Kakar  S,  Einhom  TA,  Gerstenfeld 
LC.  Genetic  Diversity  of  Skeletal  Ontogeny  Is  Recapitulated  During  Fracture  Healing.  Presented  at  the 
27th  Annual  Meeting  of  the  American  Society  for  Bone  and  Mineral  Research.  Nashville,  Tennessee 
USA.  September  23-27,  2005.  (POSTER)  -  Appendix  1. 

Stapleton  S,  Maiouf  D,  Hindson  D,  Jepsen  KJ,  PriceC,  Gerstenfeld  LC.  Genetic  Variability  in  Osteoclast 
Activity  and  3-Dimensional  Spatial  Distribution  Within  14  Day  Post  Birth  Femurs  Presented  st  the 27th 
Annual  Meeting  of  the  American  Society  for  Bone  and  Mineral  Research.  Nashville,  Tennessee,  USA 
September  23-27,  2005  (POSTER)  -  Appendix  2. 

Stapleton  S,  McLean  J,  Maiouf  D,  Hindson  D,  PriceC,  Jepsen  KJ,  Gerstenfeld  LC.  Genetic  Variation  in 
theOrteoclast  Activity  and  3>atial  DistribUion  During  Growth  and  Modeling  of  Long  Bones  as  Assessed 
by  3-Dimensional  Reconstruction.  Presented  si  the  52^  Annual  Meeting  of  the  Orthopaedic  Research 
Society.  Chicago,  Illinois,  USA.  March  19-22,  2006.  (POSTER)  -  Appendix  3. 

Jepsen  KJ,  PriceC,  Nasser  P,  Hu  B,  Alapatt  M,  Slkman  L,  Sapleton  S,  Kakar  S,  Einhom  TA, 

Gerrtenfeid  LC.  Genetic  Variability  in  the  Rate  of  Fracture  Healing  Within  Inbred  Srainsof  Mice 
Reveals  Basic  Genetic  Differences  in  the  Rates  of  Epiphyseal  Growth.  Presented  at  the  52nd  Annua] 
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Meeting  of  the  Orthopaedic  Research  Society.  Chicago,  Illinois,  US^.  March  19-22, 2006.  (POSTER)  - 
Appendix  4. 

Papers 

Wang  K,  Vishwanath  P,  Eichler  G,  Edgar  C,  Bnhom  TA,  Smith  T,  Gerstenfeld  LC.  Analysisof  Fracture 
Healing  by  Large  Scale  Transcriptional  Profile  Identified  Temporal  Relationships  Between 
M etalloproteinase  and  ADAM TSmRNA  Expression.  Matrix  Biology,  submitted.  -  Appendix  5. 

Gerdenfeld  LC,  Alkhiary  YM,  Krall  EA,  NicholIsFH,  Stapleton  SN,  Fitch  JL.Tsa/  AW,  Graves DT, 
Jepsen  KJ,  Bnhom  TA.  Three  Dimensional  Reconstruction  of  Fracture  Cal  I  us  Morphogenesis 
Demonstrate  A^mmetry  in  Callus  Development.  Journal  of  Bone  and  Mineral  Resrerrh,  submitted.  - 
Appendix  6. 


There  is  a  minimum  of  four  full  length  articles  that  will  be  written  in2006that  are  derivative  from  this 
work.  I  wi  1 1  be  happy  on  request  to  submit  these  compl  eted  manuscri  pts  and  f  i  nal  dtati  ons  of  thi  s  work 
as  it  is  published. 

Conclusions 


1.  Genetic  variability  in  bone  quality  is  recapitulated  during  bone  healing. 

2.  Genetic  variability  in  bone  quality  is  relatableto  variations  in  rates  of  bone  heating. 

3.  Genetic  variability  in  rates  of  bone  healing  appears  to  be  due  to  genetic  differences  in  endochodrna 
development  a  growth. 

4.  The  identification  of  spedfictypes  of  struct urd  and  biomechanical  traitsthat  are  both  relatefcleto  rdes 
of  fracture  healing  and  spedfic  biological  processes  that  occur  during  fracture  healing  have  the  potent  id 
in  an  intermedide  period  of  time  (2-five years)  to  be  developed  intodinicdly  usable  assays  thd  are 
prognodic  of  slow  or  fast  fracture  heding  and  of  d  risk  to  fracture  bone  qudity  to  stress  fractures 
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Appendix  1 

Genetic  Diversity  of  Skeletal  Ontogeny  Is  Recapitulated  During  Fracture  Healing 

K.  J  Jepsen1,  C.  Price1,  P.  Nasser*1,  M.  Alapatt*2,  a  N.  Stapleton2,  L.  Slkman*2,  S.  Kakar2,  T.  A. 
Enhom  ,  L.  C.  Gerstenfeld 

1  2 

Department  of  Orthopaedics,  Mount  Snai  School  of  Medicine,  New  York,  NY,  USA,  Orthopaedic 

Surgery,  Boston  University  Medical  Center,  Boston,  MA,  USA. 

Microarchitecture  (mi croCT)  and  mechanical  properties  (torque  to  fai I ure)  were  examined  over  35  da/s 
of  fracture  heal  i  ng  i  n  A/J,  C57BL/6J  (B6) ,  and  C3H/HeJ  (C3H)  i  nbred  strai  ns  of  mi  ce.  These  studi  es 
demonstrated  that  the  major  architectural  characteristics  of  bone  quality  (moment  of  inertia  cortical  area, 
trabecular  volume),  which  varied  significantly  among  the  strains,  were  recapitulated  during  fracture 
healing.  The  strain  having  the  highest  pre-fracture  strength  (C3H)  however  regained  this  property  the 
slowed.  Both  micro-CT  analysis  and  thetemporal  expression  patterns  of  cartilage  versus  bone  mRNAs 
demonstrated  very  slow  (>2wk  delay)  cartilage  resorption  for  C3H  compared  to  B6.  This  was  due  to  the 
slower  transit  of  C3H  bones  through  the  endochondral  phase.  Furthermore,  this  si  ewer  transit  time  could 
be  related  to  genetic  variability  in  the  rate  of  endochondral  maturation  based  on  the  5-7  day  difference  in 
times  of  peak  chondrocyte  hypertrophy  based  on  Col  10A1  expression.  Quantitative  analysis  of  the 
expression  of  both  cartilage  and  bone  extracellular  matrix  mRNAs  showed  a  direct  correlation  between 
thetotal  anabolic  activity  of  a  given  strain  of  mouse  and  the  overall  volume  of  bone  tissue  formed  such 
that  (C3H>B6>A4).  Examination  of  14  day  post  birth  growth  plates  demonstrated  similar  temporal 
patterns  of  endochondral  progression  with  the  AJ  and  B6  mouse  strains  showing  substantially  advanced 
stages  of  maturation  as  compared  to  C3H  that  lagged  behind.  These  data  suggested  that  the  genetic 
elements  which  control  the  primary  structural  attributes  of  bone  quality,  are  recapitulated  during  fracture 
healing  and  can  be  easily  assessed  in  thetemporal  window  of  fracture  healing.  These  data  indicate  that  the 
basic  temporal  and  quantitative  attributes  of  endochondral  bone  formation  during  fracture  healing  ae 
genetically  regulated.  Finally,  we  show  that  the  rate  of  bone  healing  is  not  simply  related  to  varied  I  ity  in 
bone  quality  but  is  also  controlled  by  other  factors  associated  with  genetic  variability  intrinsic  to 
endochondral  bone  formation. 
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Appendix  2 


Genetic  Variability  in  Osteoclast  Activity  and  3-Dimensional  Spatial  Distribution  Within  14  Day 
Post  Birth  Femurs 

S.  Sapleton*,  D.  Malouf,  D.  Hindson,  K.  J.  Jepsen,+  C.  Price*  LC Gerstenfeld* 

'Orthopaedic  Research  Laboratory,  Boston  University  Medical  Center,  Boston,  MA  02118-2526,  USA, 
Department  of  Orthopaedics,  Mount  Snai  School  of  Medicine,  New  York,  New  York  10029-6574,  USA. 

Numerous  studies  have  shewn  that  variation  in  bone  morphology  in  different  inbred  strains  of  mice  is 
genetically  determined.  The  genetic  variation  is  reflected  in  different  rates  of  formation  aid  resorption  on 
endorteal  vs.  periosteal  surfaces;  the  relative  amount  and  location  of  cellular  activities  thus  determining 
the  relationship  between  bone  size  and  the  mechanical  demands  imposed  by  weight  bearing  during 
growth.  A  quantitative  measure  of  osteoclast  numbers  and  their  3-dimensional  spatial  distribution  were 
assessed  in  A/J,  C57BL/6J  (B6),  and  C3H/HeJ(C3H)  strains  of  mice  Transverse  serial  sections  were 
taken  at  100  micron  increments  along  the  length  of  femurs  from  14  day  mice  and  stained  for  osteoclast 
activity.  Serial  sections  were  arranged  in  stacks  and  the  3-D  spatial  distribution  of  osteoclasts  within  the 
femurs  was  reconstructed  (A  J  and  C3H  are  presented) .  Renderi  ngs  of  the  bone  i  n  differi  ng  I  ongitudi  nd 
planes,  or  for  TRAP  activity  alone,  were  examined.  Moving  distaily  along  the  femur  revested  an 
approximate  50/50  distribution  of  TRAP  staining  on  both  endosteal  and  periosteal  surfaces  within  A  J  aid 
B6  femurs,  but  no  discernable  pattern  for  C3H.  The  AJ  strain,  which  has  the  smatlest  bone  mass,  had  the 
higher  level  of  TRAP  activity  with  the  majority  of  the  enzyme  activity  seen  either  on  the  periosteal  or 
opposing  endorteai  surfaces  (compae  the  distribution  of  the  pattern  of  dak  grey  coloring  that  denotes 
aeas  of  TRAP  activity).  For  the  B6  strain,  which  had  the  next  highest  level,  most  of  the  activity  was 
ether  intracortical  or  on  trabecula  bone  in  the  medullary  space;  much  I  ess  staining  visualized  on  the  bone 
surfaces  Intaertingly,  this  strain  has  the  thinnest  cortices  and  I  agest  medullary  space.  TRAP  staining  in 
theC3H  animals  was  almost  unobservable  when  viewed  in  the  reconstructions  and  this  strain  has  both  the 
thickest  cortices  and  highest  bone  mass.  These  data  demonstrate  profound  differences,  in  both  the 
quantity  aid  spatial  chaacteri sties  of  osteoclastic  activities,  in  different  strains  of  inbred  mice.  This 
aiggerted  that  both  these  features  of  cel  I  ul  a  acti  vity  ae  rel  atabl  e  to  spati  al  geometry  and  overd  I  mass  of 
a  bone  and  are  genetically  regulated. 
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Appendix  3 


GENETIC  VARIATION  IN  THE  OSTEOCLAST  ACTIVITY  AND  SPATIAL  DISTRIBUTION  DURING  GROWTH  AND  MODELING  OF 
LONG  BONES  AS  ASSESSED  BY  3-DIMENSIONAL  RECONSTRUCTION 

♦‘Stapleton,  S;  ‘McLean,  J;  ‘Malouf,  D;  ‘Hindson, D;  “Rice,  C;  “Jepsen,  KJ;  ‘Gerderfeld,  LC 
‘Orthopaedic  Research  Laboratory,  Boston  University  Medical  Center,  Boston,  MA  02118-2526,  USA 
“Department  of  Orthopaedics,  Mount  Snai  School  of  Medicine,  New  York,  New  York  10029-6574,  USA. 

aistaple@bu.edu 


INTRODUCTION 
Bone  quality  has  been  shown  to  be  control  led  by  multiple  genes  that 
regulate  skeletal  tissue  structure.  One  set  of  genes  dfects  geometri cat 
features  such  as  cortical  and  trabecular  thickness,  cross  sectional  areas 
and  trabecular  numbers  whilea  second,  not  necessity  independent,  set 
of  genes  affects  tissue  material  properties  inclusive  of  mineral  content 
and  organic  matrix  composition  (1).  At  the  cellular  level  these  genetic 
variations  are  reflected  in  different  rates  of  formation  and  resorption  on 
endosteal  and  periosteal  surfaces  versus  amount  of  remodeling  in 
oorticd  and  trabecular  bone.  Thus,  the  relative  amount  and  location  of 
cellular  activities  isa  large  determinant  of  many  of  the  geometric 
features  of  skeletal  tissues  and  is  responsive  to  the  mechanical  demands 
i  mposed  by  weight  bearing  duri ng  growth.  I  n  this  study  a  quantitative 
measure  of  the  spatial  distribiiion  of  osteodast  activity  was  assessed  in 
three  strains  of  inbred  male  mice  (A/J,  C57BL/6J  [B6],  and  C3H/HeJ 
[C3H]), which  have  varying  bone  qualities.  The  distribution  of 
osteoclast  activity  was  oompared  relative  to  the  overdl  endochondral 
and  osseous  growth  of  these  ti  ssues  as  detemi  ned  by  three  dimensi  oral 
histological  reconstruct  on  of  the  femurs  of  these  mi  ce. 

MATERIALS  AND  METHODS 

R ve  micron  transverse  secti ons  were  generated  from  i ntact  fourteen  day 
decalcified  and  paraffin  embedded  mouse  femurs.  Sections  were  taken  at 
100  micron  increments  across  the  entire  length  of  the  femur  aid 
saf  ranin-O  fast  green  staining  was  used  to  dstinguish  mature  cartilage 
from  boneand  other  connective  tissues.  Osteodast  activity  was  detected 
by  the  use  of  TRAP  staining.  Sections  were  then  photographed  at  40X 
and  loaded  into  Adobe  Fhotoshop  for  editing.  1  mages  were  aligned  with 
each  other  and  pertinent  tissues  identified.  Boneand  cartilage  were 
di  dingui  shed  by  hand  tracing,  and  TRAPstai  ned  cells  were  located  by 
pixel  value.  All  sections  were  false  colored  according  to  type  of  tissue 
that  was  being  assessed  aid  then  converted  into  grey  scale  Amira  3.0 
software  system  (Visual  Concepts  GmbH ,  Company,  Konrad-Zuse- 
Zentrum(ZIB),  Research  Institute  Bert  in,  Germany)  was  used  for  3D 
visiBliz^ion  and  data  analysis.  All  image  modifications  were 
performed  in  Adobe  Photoshop  Software,  (Adobe Systems, 

Incorporated;  San  Jose,  CA  USA).  Areas  of  TRAP  activity  were 
calculated  ty  use  of  the  surface  area  function  in  Amira  3.1  within  each 
of  the  reconstructed  models.  Areas  and  spatial  distribution  of  TRAP 
activity  was  further  analyzed  by  selective  cell  counts  of  individual 
osteoclasts  on  specific  anatomical  surfaces. 

RESULTS 

Analysisof  the  three  strains  of  femurs  showed  that  there  was  an 
approximate 50/50  distribution  of  TRAPstaining  on  both  endosteal  and 
periosteal  surfaces  within  AJ  and  B6  femurs  within  the  center  di^ihysaal 
regions  of  the  bone,  but  that  no  discemable  pattern  was  seen  for  C3H 
(Rgurel).  TheAJstrain,  which  has  the  smallest  bone  mess,  had  the 
highed  level  of  TRAP  activity  with  the  majority  of  the  enzyme  activity 
seen  either  on  the  periosteal  or  opposing  endosteal  surfaces  (corrp»e  the 
distribution  of  the  pattern  of  dark  grey  coloring  thd  denotes  areas  of 
TRAPactivity).  For  the  B6  strain,  which  had  the  next  highest  level,  a 
pattern  amilarto  AJ  resulted,  but  with  less  staining  visualized  on  the 
bone  surfaces  Interestingly,  thisstrain  has  the  thinnest  oorticesand 
larged  medullary  space  with  the  greeted  number  of  cells  seen  in  the 
intramedullary  spaoe.  TRAP  staining  in  the  C3H  animalswas  aimed 
unobservable  when  viewed  in  the  reconstructions  and  this  strain  has  both 


the  thick ed  cortices  and  highest  bone  mass.  An  analysis  of  the  actual 
TRAP  da  ned  area  (Figure  2)  and  cells  also  showed  varying  quantitative 
results  between  the  three  drains  of  mice 

An  examination  of  the  progression  of  endochondral  tissue 
maturation  in  the  three  drains  demondrated  that  the  AJ  stains  had  the 
mod  progressed  endochondral  mduration,  B6  strain  was  intermediate 
and  the  C3  had  the  lead.  In  situ  hybridization  with  Type  II  oollagen 
validated  strain  specific  variation  in  chondrocyte  metabolic  activity,  and 
an  assessment  of  cellular  numbers  within  growth  plate  col  urms  showed 
that  the  A  J  strai  n  had  the  greatest  cel  I  number  per  column  aid  larged 
column  heights  while  the  C3  mice  had  the  smal  led  number  of  cells  and 
shorted  columns. 

DISCUSSION 

These  dda  demonstrate  profound  differences,  in  both  the  quantity  and 
spatial  charaderidics  of  odeodastic  activities,  in  different  strains  of 
inbred  mice.  These  data  sugged  that  both  these  features  of  odeodad 
activity  are  relatable  to  spatial  geometry  and  overall  mass  of  a  bone  and 
are  genetically  regulated.  Data  examining  the  epiphyseal  growth  plates 
suggeded  that  there  isalso  genetic  variation  in  the  rate  of  linear  growth 
that  isassodated  with  the  rate  of  endochondral  chondrocyte  maturdi on 
such  that  the  AJ  showed  the  faded  rate  of  growth  and  cartilage 
maturation  while  the  C3H  had  the  slowed  rate  and  cartilage  mstwation. 
These  dda  pose  speafic  quedions  as  to  whether  there  is  an  assoaation 
between  the  smal  ler  transrerse  bone  areas  and  the  rate  of  I  iner  bone 
growth,  and  what  are  the  molecular  mechanisms  of  linkage  between 
transverse  and  linear  growth. 
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Fig  1.  Spatial  Orientation  of  modeling  and  endochondral  cartilage 
development  in  inbred  drains  of  mice.  Three  views  (left  to  right)  of  the 
recondruded  spdial  didribution  of  TRAP  (odeodast)  activity  in  the 
bone  in  the  3  strai  ns  of  mice  The  bone  is  modeled  in  light  gray  aid  the 
TRAPactivity  is  modeled  in  black.  The  left  view  depicts  the  periosteal 
surface.  The  middle  view  is  a  cut  away  looking  onto  theendodeal 
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surface.  The  right  view  depicts  the  enti  re  volume  that  was  sampl  ed  for 
TRAP  staining  without  the  superposition  of  the  reconstruction  of  the 
bone  surfaces,  The  far  right  panel  modelsthe  mature  cartilage  (light 
grey)  and  immature  cartilage  (blade)  within  the  growth  plates 
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Fig  2.  Catabolicactivity  per  bonaareainAJ,  B6andC3H  mice 
Area  calculations  based  on  TRAP  staining  in  three  strains  of  mice. 
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INTRODUCTION 

Fracture  healing  isaspedalized  post-natal  repair  process  that 
recapitulates  many  aspects  of  embryological  skeletal  development. 
Recent  studies  using  inbred  strains  of  mice  have  shown  that  bone  quality 
is  controlled  by  multiple  genes  that  regulate  skeletal  tissue  structure. 
Because  many  of  the  molecular  mechanisms  that  control  cellular 
differentidion  and  growth  during  embryogenesisare  recapitulated 
during  fracture  repair,  we  tested  the  hypothesis  that  the  rate  of  fracture 
heeling  will  vary  with  genetic  background.  Wefurlher  tested  whether 
genetic  variability  infracture  repair  may  have  its  on  gins  in  variable 
endochondral  boneformetion  processes  by  comparing  the  cellular 
behavior  during  fracture  repair  with  the  cellular  behavior  of  the  growth 
plate.  In  this  study  three  drainsof  inbred  mice  with  varying  bone 
qualities  (A/J,  C57BL/6J[B6],  and  C3H/HeJ  [C3H])  were  fractured  and 
the  rate  of  regain  of  mechanical  properties  was  assessed.  Basic 
histological  and  moleaiar  features  across  fracture  heali  ng  were  used  to 
measure  rates  of  chondrogenic  aid  osteogenic  differentiation.  These 
studies  reveeded  that  the  overall  anabolic  activity  of  skeletal  tissues  is 
genetically  controlled  and  that  rates  of  chondrocyte  maturation  during 
endochondral  bone  formation  also  show  genetic  variation  that  affect 
rates  of  healing.  These  observe ons  led  to  comparisons  of  the  rates  of 
chondrogenic  differentiation  during  post-natal  epiphyseal  bone  growth 
that  validated  basic  genetic  differences  in  skeletal  tissue  development. 
MATERIALS  AND  METHODS 

Srrple  transverse  mid-diaphyseal  femur  fractures  were  generated  in  8  to 
10  week  old  male  A/J,  C57BL/6J  [B6] ,  and  C3H/HeJ  [C3H]  trace  as 
previously  described  (1).  Fracturecailuses  were  harvested  at  21  raid  35 
daysfor  torsion  mechanical  testing  and  microCT  analysis.  Callus  tissues 
were  harvested  at  multiple  post-fracture  time  points  over  a  28  day  period 
for  mRNA  and  histological  measurements.  The  respective  distal  and 
proximal  epiphyseal  growth  platesof  tibia andfemurs from  14 day  old 
(pod  birth)  mce  were  histologically  examined.  Type  II  and  X  collagen 
mRNA  expression  analysis  was  assessed  by  in  situ  hybridization  and 
mRNA  expression  within  the  whole  tissue  was  assessed  by  ribonudease 
protection  analysis 
RESULTS 

Microarchitecture  (microCT)  and  mechanical  properties  (torqueto 
failure)  were  examined  at  21  and  35  days  of  fracture  healing  in  A/J, 
C57BL/6J(B6),  and  C3H/HeJ  (C3H).  These  studies  demonstrated  that 
the  major  architectural  characteristics  of  bone  quality  (moment  of 
inertia,  cortical  area,  trabecular  volume),  which  varied  significantly 
among  the  strains,  were  recapitulated  during  fracture  healing  (Figure  1). 
The  strain  having  the  highest  pre-fracture  strength  (C3H)  however 
regained  its  strength  the  slowest  as  based  on  torque  to  failure 
measurements.  Both  rra'cro-CT  analysis  (seen  as  x-ray  lucency)  and  the 
temporal  expression  patterns  of  cartilage  versus  bone  mRNAs 
demonstrated  very  slow  £  2  wk  delay)  cartilage  resorption  for  C3H 
compared  to  either  B6  or  A  J.  This  was  due  to  the  slower  transit  of  C3H 
bones  through  the  endochondral  phase  based  on  the  ratios  of  Col2a1  to 
CollOal  expression  that  were  used  as  a  measure  of  the  progresion  of 


hypertrophic  chondrocyte  differentiation  (Figure  2).  Furthermore,  this 
slower  transit  time  could  be  related  to  genetic  variablity  in  the  rate  of 
endochondral  msturati on  based  on  the  5-7  day  difference  in  times  of 
peak  chondrocyte  hypertophy  based  on  CollOAl  expression. 
Quantitative  analysis  of  the  expression  of  both  cartilage  and  bone 
extracellular  matrix  mRNAsfurther  showed  a  drect  correlation  between 
the  total  anabolic  activity  of  a  given  strain  of  mouse  and  the  overall 
volume  of  bone  tissue  formed  such  that  (C3H>B6>AJ>. 

Examination  of  14  day  post  birth  growth  plies  demonstrated  similar 
temporal  patterns  of  endochondral  progression  with  the  A  J  and  B6 
mouse  strains  showing  substantially  advanced  stages  of  nratiriion  as 
compered  to  C3H  which  lagged  behind.  This  is  seen  in  Figure  3  of  the  in 
situ  hybridization  of  the  femoral  growth  plates  that  showed  the 
associated  genetic  variablity  in  overall  Col2a1  expression  intrinsic  to  the 
deferring  strains  c#  mice.  The  progression  of  endochondral  tissue  was 
markedl  y  different  i  n  that  A  J  stai  ns  had  the  most  progressed 
endochondral  maturation  with  B6  strain  bang  intermediate  and  the  C3 
had  tiie  least.  Cellular  numbers  within  growth  plate  oolumnsalso 
showed  that  the  AJ  strain  had  the  greatest  cell  number  per  column  and 
largest  column  heights  while  the  C3  mice  had  the  smallest  number  of 
cells  and  shortest  columns, 

DISCUSSION 

These  data  suggested  that  the  genetic  el  ements  which  control  the 
primary  structural  attributes  of  bone  quality  are  recapitulated  during 
fracture  healing  and  can  be  easily  assessed  in  the  temporal  window  of 
fracture  healing.  These  data  indicate  that  the  basictemporat  and 
quantitative  attributes  of  endochondral  bone  formation  during  fracture 
healing  are  genetically  regulated.  Finally,  we  show  that  the  rate  of  bone 
healing  isnot  simply  related  to  variability  in  bone  quality  butisalso 
controlled  by  other  factors 
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Fig  I.  Structural  and  Biomechanical  Comparisons  of  Fracture  Cell  uses 
From  Three  Inbred  Strains  of  Mice.  Representative  microCT  imagesin 
both  transverse  and  longitudinal  orientations  taken  at  21  and  35  days 
post  fracture.  The  orientation  of  the  microCT  images  is  denoted  for  the 
I  aft  most  transverse  section  of  the  A  J  strain.  Boxed  area  of  C3H  day  21 
image  shows  area  that  remained  unbridged  by  bone. 
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Fig  2  Ratios  of  Bone  to  Cartilage  and  Hypertrophic  Cartilage 
Maturiion  Over  Time  Course  of  Fracture  Cai  I  us  Healing. 

The  measured  ratios  of  incS  virtual  collagen  mRNA  expression  are 
denoted inthefigure.  Solid line=C3H,  Dotted  line AJ,  DashedlineB6. 
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Fig  3  Comparison  of  Distal  Fenxir  Growth  Hates  at  14  Days  Post  Birth. 
Upper  panel  are  li^tt  level  images  of  representative  areas  from  the 
growth  platesof  thethreestrainsof  mice.  Lower  panels  showed  dark 
field  images  from  the  in  situ  hybridzations  with  Col2A1  probes, 
Staggered  arrows  show  the  variation  i  n  col  urm  sizes  of  the  hypertrophic 
and  growth  zones  of  the  three  strains.  Numbers  are  simple  counts  of 
cells  seen  in  the  two  zones 

1.  Kon  etal  .  J  Bone  Miner  Res  2001;  16:1004-1014. 
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Appendix  5 


Analysis  of  Fracture  Healing  by  Large  Scale  Transcriptional  Profile  Identified  Temporal  Relationships  Between 

M  etal  I  oprotei  nase  and  A  DAM  TS  mRNA  Express  on 
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Abstract 


The  a  m  of  thi  s  study  was  to  val  i  date  the  use  of  transcri  pti  onai  prof  i  I  i  ng  as  a  means  of  characteri  zi  ng  the  compi  ex 
i  nteracti  ons  of  the  thousands  of  genes  that  are  expressed  duri  ng  fracture  heal  i  ng.  Standard  mi  d-diaphy9eal  ti  bi  a 
fractures  were  generated  in  C57/B6  murine  tibiae  and  the  transcriptional  expression  of  -13,000  genes  was 
assessed.  Three  time  points  after  fracture  were  assessed:  d^  three,  representative  of  the  inflammatory  phase;  day 
10,  representdive  of  the  peak  of  cartilage  formation;  and  d^  21,  representative  of  the  period  of  primary  bone 
formation  and  coupled  remodeling.  A  self-organizing  mapping  approach  of  the  data  revealed  the  temporal 
relationships  between  the  expression  of  mRNAs  for  extracellular  matrix  proteins  and  the  proteases  that  degrade 
the  proteoglycan  and  collagenous  matrices  A  broad  group  of  extracellular  matrix  protein  mRNAs  representative 
of  basement  membranes,  blood  vessels  ard  cartilage  all  showed  elevated  expression  over  the  first  21  d^sof 
fracture  healing.  The  sorting  of  the  data  identified  an  orderly  temporal  expression  of  the  metal  I  oproteinases  and 
ADAM TS  during  the  progressi  on  of  fracture  heal i ng  with  (M  M  P2/M  M  P14/TI M  P2)  and  ADAM TS  4  and  1 5 
preceding  the  expression  of  (MMP9/MMP13).  Based  on  their  patterns  of  expression,  relative  to  the  known 
activities  of  the  encoded  proteolytic  enzymes,  our  results  suggest  that  the  dissolution  of  cartilage  protoeglycans 
proceeds  before  the  underlyi  ng  col  I  agenous  components  of  the  mdrix  are  removed.  The  exd  ua  on  of  several 
mRNAs  that  are  normally  expressed  by  osteoclasts  in  the  profiles  of  mRNAs  from  days  3  and  10  suggests  that 
osteoclastic  activity  was  largely  absent  during  the  early  periods  of  cartilage  tissue  formation  andthd  proteoglycan 
and  specific  collagenase  activities,  precedes  or  is  prerequisite  to  later  osteoclast  infiltration  into  the  remodeling 
tissues. 

Keywords:  MMP,  ADAMTS,  Fracture  Healing,  Transcriptional  Profiling 
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1.  Introduction 


Osseous  skel  etal  ti  ssues  are  one  of  the  few  ti  ssues  i  n  adult  vertebrate  ani  mals  that  regenerate  i  n  response 
to  injury.  Fracture  healing  encompasses  the  complex  interaction  of  multiplecell  types  and  many  different 
biological  processes  (Ferguson  et  al„  1999;  Gerstenfeid  et  al.,  2003b;  Vortkamp  et  al.,  1998).  Thus  fracture 
heal  i  ng  i  s  a  uni  que  model  to  i  dentify  the  rel  ati  onshi  ps  between  the  mechani  sms  thd  regul  ate  sped  fi  c  responses  to 
atraumatic  injury  and  the  genetic  components  that  control  the  development  of  bone.  The  identification  and 
understanding  of  interactions  of  the  various  genomic  components  th^  both  respond  to  injury,  affect  bone 
development  and  regul  ate  repair  are  therefore  of  considerable  dini  cal  importance.  Identifying  these  molecular 
mechanisms  help  us  develop  predictive  indices  for  the  rates  of  fracture  healing  that  are  based  on  spedfic 
biological  processes.  These  data  will  also  be  informdtiveto  potential  therapeutic  interventions  that  may  promote 
bone  healing. 

M  any  previ  ous  studi  es  of  fracture  heal  i  ng  have  focused  on  the  activiti  es  of  i  ndi vi  dual  genes  that  are  of 
functional  importance  during  the  differentiation  of  skeletal  cell  types,  The  biological  complexity  of  fracture 
healing,  however  suggests  that  it  is  regul  ded  by  hundreds,  if  not  thousands;  of  transcriptionally  regulated  genes. 
These  genes  are  expressed  by  both  skeletal  and  non  skeletal  cell  types  and  they  are  involved  in  multiple  biological 
processes  that  fad  I  itate  bone  healing.  In  order  to  understand  how  the  adivities  of  such  a  large  number  of  genes 
are  coordinated,  and  how  their  activities  are  assodated  with  a  given  cellular  or  biological  process,  recent 
innovations  in  the  ability  to  analyze  genome  wide  transcriptional  activities  can  be  used  (Hadji  argyrou  et  al.,  2002; 
Desai  et  al.,  2003).  Such  approaches  enable  the  transcriptional  activities  of  large  numbers  of  genes  to  be  assessed 
and  the  complex  interactions  of  genes  to  be  organized  and  visualized  in  an  understandable  and  interpretable 
manner. 

The  ai  m  of  thi  s  study  was  to  val  i  date  the  use  of  transcri  pti  onal  prof  i  I  i  ng  as  a  means  of  characterizi  ng  the 
complex  interactions  of  approximately  13,000  genes  that  are  expressed  during  fracture  healing.  A  simplified 
study  desi  gn  usi  ng  a  I  i  mi  ted  set  of  ti  me  poi  nts  that  woul  d  del  i  neate  the  i  nfl  ammatory ,  endochondral  (carti  I  age) , 
and  primary  bone  formation  phases  of  fracture  healing  was  carried  out.  This  study  selectively  focused  on  the 
complement  of  ©dracellular  matrix  proteins  and  the  two  primary  proteolyti  c  enzyme  f  ami  I  i  es  (metal  I  oprotei  nases 
and  ADAMTs)  that  are  responsible  for  remodeling  extracellular  matrices.  This  analysis  identified  several 
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hundred  expressed  genes  that  showed  significant  differential  expression  compared  to  unfractured  bones. 
Clustering  analysis  identified  specific  groups  of  metalloproteinases  and  ADAMTS  mRNAsand  placed  their 
expression  in  context  to  the  expression  of  specific  mRN As  associated  with  cartilage,  bone  and  osteoclast  cell 
functions. 

2.  Results 

The  expression  of  a  group  of  extracellular  matrix  protdn  mRNAsthat  is  representative  of  cartilage  and 
bone  formation  was  first  assessed  in  each  set  of  RN  As  that  was  used  for  transcriptional  profiling.  Ribonudease 
protecti  on  analysis  was  carri  ed  out  i  n  order  to  ensure  that  each  RN  A  sampl  e  was  representati  ve  of  the 
inflammatory,  cartilage  and  bone  formation  phases  of  fracture  healing  (Figure  1).  The  inflammatory  phase  of 
fracture  heal  i  ng  (d^  3) ,  was  characteri  zed  by  an  i  niti  al  decrease  i  n  the  expressi  on  of  the  maj  or  mRN  A  markers  of 
differentiated  osteogenic  ceils  (osteocalcin  type  I  collagen  and  bone  a  a  oprotein).  As  previously  shown,  this  is 
due  to  the  large  influx  of  non  osteogenic  cells  within  the  hematoma  that  forms  in  response  to  the  fracture  (Kon  et 
a.,  2001).  At  thistime  large  numbers  of  skdetogenic  precursors  tha  contribute  to  the  i  niti  a  repar  tissue  are 
recruited  but  do  not  yd  express  a  skdeta  tissue  phenotype.  Consistent  with  this  interpretdion  was  the  increased 
expression  of  osteopontinthd  is  expressed  by  macrophages  and  other  inflammdory  cells  (Kon  et  a.,  2001),  as 
wdl  as  early  skdeta  genic  precursors  (Barnes  et  a.,  1999;  Bolander,  1992).  Da/  10  ca  I  uses  tissue  represents  the 
peak  of  cartilage  formaion  thd  takes  place  (tiring  endochondra  boneformaion(Gerdenfdd  da.,  2003a).  The 
redricted  ndure  of  the  chondrocyte  mRNA  expression  d  this  time  point  is  exemplified  by  the  high  levds  of  both 
collagen  types  1 1  andX  mRNA  expression.  D^21  is  representdive  both  of  primary  bone  formdion  and  coupled 
remodding.  At  thistime  point,  the  profile  of  expressed  mRN  As  has  d  most  returned  to  thd  seen  in  unfractured 
bone  however  the  levds  of  osteocddn  remdn  considerably  devded  above  the  basdine  of  normd  bone. 

Whilethe  RPA  andysis  is  an  important  assessment  of  RNA  integrity  and  provides  a  means  of  confirming 
thetempord  progression  of  the  repdr  process,  microarr^  experiments  necessitde  a  much  higher  purity  of  RNA 
thd  may  only  be  empirically  ddermined  by  carrying  out  an  actud  microarray  experiment.  The  raw  data  from  the 
spot  intensities  for  the  expression  of  osteocdcin  was  chosen  to  vdidde  the  integrity  of  the  profiling  experiments 
The  expression  of  this  mRNA  is  known  to  be  one  of  the  most  highly  restricted  to  bone,  and  its  expressi  on  has 
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been  shown  to  progressively  increase  over  the  time  course  of  fracture  hewing.  The  spot  intensities  for  both 
fluorescent  channels  from  the  reciprocal  labeling  and  hybridization  of  one  set  of  experimental  and  reference 
samples  are  presented.  These  data  show  an  excellent  quantitative  correspondence  of  the  increasing  intensity  of 
fl  uorescent  readout  over  the  ti  me  course  of  fracture  heal  i  ng,  i  n  compari  son  to  the  RPA  assay. 

Using  unfractured  bone  mRNA  as  the  reference  to  which  experimental  samples  were  compared,  we  used 
a  conservative  selection  process  in  which  only  those  mRNAs  showing  a  >log2  difference  in  expression  compared 
tothelevelsof  expression  in  unfractured  bone  would  be  considered  for  further  analyse  The  global  representation 
of  these  data  is  summarized  as  a  ratio  intensity  plot  that  is  depicted  in  Figure  3.  This  figure  illustratesthe  total 
number  of  observed  mRNAs  at  each  time  point  having  a  >log2  fold  variation  from  the  observed  leveisin 
unfractured  bone.  The  conformance  of  the  majority  of  the  spots  within  1  SO  of  their  log  ratios  a  so  provides  an 
excellent  index  of  the  overall  integrity  of  the  numeric  data  that  was  developed  from  each  arr^  (Tamayo  et  al., 

1 999;  Kohonen  and  Somervuo,  2002) .  A  compari  son  of  the  three  graphs  shows  that  the  greatest  number  of  genes 
having  >log2  variation  lesser  or  greater  expression  from  unfractured  bone  were  seen  at  days  3  and  10  d^s  post 
fracture  Far  fever  genes  showed  such  I  argequantiteive  differences  in  their  expression  relative  to  urrfractured 
bone  by  21  days  post  fracture  It  is  interesting  to  note  that  asthe  fracture  healed,  the  complement  of  genesthat 
were  differentially  expressed  greatly  diminished  asthe  bone  returnsto  its  pre  fracture  makeup  (Figure  3).  Thus 
the  complexity  of  the  number  of  biological  processes  that  are  concurrently  operative  during  the  early  phases  of 
fracture  healing  is  much  greater  than  that  which  occurred  during  the  later  periods  of  primary  bone  formation 
The  first  level  of  analysis  was  to  identify  those  expressed  mRNAs  that  are  quantitatively  different  from 
unfractured  bone  A  small  subset  of  all  of  the  mRNAs  showing  differential  levelsof  expression,  are  presented  in 
Table2.  The  primary  focus  of  thistable  isto  highlight  expressed  mRNAsthd  are  selectively  expressed  during  the 
three  phases  of  fracture  healing,  and  that  are  known  to  be  associated  with  either  cartilage  formation  or  osteoclast 
activity.  A  number  of  ether  mRNAsthat  encode  specific  morphogenetic  regulators  and  transcriptional  factors  are 
at  so  hi  flighted  in  the  tdde.  These  were  included  si  nee  they  have  not  been  previously  well  characterized  during 
fracture  healing  but  showed  strong  differences  in  thei  r  expression  relative  to  unfractured  bone  and  may  be  of 
general  interest. 
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At  day  10  post  fracture,  very  large  changes  in  the  expression  of  perlecan,  cartilage  oligomeric  matrix 
protein  (COMP)  and  matrix  glaprotein  (MGP)  mRNAS  were  observed.  All  three  of  these  mRNAs  or  their 
proteins  are  known  to  be  highly  expressed  in  developing  cartilage  but  not  in  normal  bone  Their  identification  at 
thistime  point  is  consistent  with  the  RPA  data  for  collagen  Types  II  and  X  and  further  validates  the  restricted 
expression  of  cartilage  associated  mRNAs  at  thistime  point.  On  the  other  hand,  a  number  of  enzymes  such  as 
carbonic  anhydrase2  and  tartrate  resistant  add  phosphatase  that  are  all  known  to  be  expressed  at  high  levels  by 
oSeodasts  (Roussel I e  and  Heymann,  2002),  were  significantly  diminished  from  their  reference  values  in  native 
bone  at  da/s  3  and  10  post  fracture.  Thediminished  expression  of  marker  mRNAsforosteodasts  was  also 
accompanied  by  the  eieviion  of  two  separate  antagonist  to  RANK-L  (major  regulator  of  osteoclast 
differentiation),  osteoprotegerin,  and  TNF  super  family  member  22.  A  number  other  mRNAs  of  interest  that  were 
induced  early  and  persisted  through  out  the  time  course  of  fracture  healing  indude  osteoblast  spedficfador  2 
(fasddin  l-like)  OSF-2  and  osteoglydn.  In  references  to  the  expression  of  OSF-2  and  osteoglydn  both  these 
genes  have  all  been  identified  in  two  separate  large  scale  in  vitro  transcriptional  profiling  studies  that  were 
focused  on  identifying  genes  within  mesenchymal  stem  cellsthat  are  regulated  by  BMPs(Balint  et  at.,  2003; 

Stock  et  al„  2004). 

The  first  stage  of  the  analysis  only  identified  candidate  genes  of  interest  based  on  changes  in  their 
quantitative  levels  of  expression.  In  the  second  stage  of  this  analysis  bioinformatics  approaches  were  ^plied  to 
identify  functional  groupings  of  genes.  A  self-organizing  mapping  (SOMs)  approach  was  used  to  identify  the 
temporal  patterns  of  expression  for  two  groups  of  proteolytic  enzymes  (M  M  PS,  ADAMTS)  and  the  extracellular 
matrix  protein  mRNAs.  Table  3  presents  the  two  groupings  of  extracellular  matrix  prots'n  mRNAsthat  showed 
slightly  different  temporal  patterns  of  expression  but  which  were  grouped  dosely  together  in  the  two  dimensional 
graphical  anaiysisof  the90Ms(data  not  shown).  Both  groups  of  these  mRNAs  show  increased  levels  of 
expression  over  unfractured  bone  at  all  of  the  time  points  The  Group  1  mRNAsshowed  an  increasing  level  of 
expression  between  3  and  14  fcys  after  which  their  levels  began  to  fall.  This  group  of  mRNAs  inducted  those 
encoding  proteins  found  mainly  in  basement  membranes  and  vascular  elements,  a  second  set  of  mRNAs  primarily 
encoding  proteins  found  in  cartilage,  and  athird  group  of  mRNAs  for  proteins  found  in  many  Afferent  connective 
tissues  Group  1 1  mRNAs  had  a  pattern  of  expression  that  showed  a  continuous  increasing  level  of  expression 
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across  the  entire  time  course  and  inducted  mRNAs  encoding  both  al  aid  a2  polypeptides  of  type!  collagen. 
Interestingly,  the  mRNAs  encoding  several  different  metalloproteinase  were  one  of  a  several  functional  groups  of 
expressed  mRNAsthat  were  identified  without  prior  selection  and  based  solely  on  their  strong  quantitative 
differences  in  expression.  The  identification  of  multiple  MM Ps  that  areinvolved  in  both  vascularization  and 
collagenous  tissue  proteolysis  was  of  considerable  interest,  s'ncethe  resorption  of  cartilage  tissue  and  tissue 
vascularization  representstwo  of  the  major  transitional  processes  of  endochondral  boneformation(Vuetal., 
1998;  Zhoetal.,  2000;  Street  etal.,  2002;  Ortega  etal.,  2003).  These  observations  led  us  to  subsequently  use 
SOM  s  to  further  exami  net  he  temporal  rel  at  i  onshi  ps  between  the  metal  I  oprotei  nases  and  the  ADAM  TS  proteases 
that  are  known  to  degrade  proteoglycans.  A  graphical  presentation  for  MM  Ps  and  TIMPs  are  presented  in  Figure 
4.  Individual  groupsof  MM P  mRNAs  showed  varying  temporal  patterns  of  expression,  with  the  expression  of 
MMPs  14  and  2  preceding  that  of  MMP9andMMP13.  In  contrast,  the  expression  of  MMPs  8  and  19  showed  a 
strong  i  nhi  biti  on  in  thei  r  expressi  on  over  the  21  days  of  bone  heal  i  ng.  The  metal  I  oprotei  nase  i  nhi  bitors  (Tl  M  Ps  2 
and  3)  had  an  almost  identical  pattern  of  expression  as  MMPs  2  and  14  while  TIMP1  was  induced  early  after 
fracture  and  remained  highly  expressed  thereafter. 

Specific  members  of  the  ADAMTS  proteases  genefamily  also  showed  unique  temporal  patterns  of 
expression  across  the  time  course  of  fracture  (Figure5).  ADAMTS  4  and  15,  which  both  have  aggrecanases 
activity  (Apte;  2004),  paralleled  the  expression  of  MMPS 2  and  14.  On  the  other  hand,  ADAMTS  1  which  also 
has  potential  aggrencanase  activity  followed  the  pattern  of  expression  of  MMPS  9  and  13.  The  final  ADAMTS 
that  we  observed  to  be  expressed  in  these  tissues  that  fit  within  the  quantitative  selection  of  log  2  or  greater 
differencewas  ADAMTS  10  which  followed  the  pattern  of  MMPs8ar»d  19. 

I  n  order  to  validiethe  mRNA  expression  patterns  of  the  M  MPs  and  ADAMTSthi  were  observed  with 
the  transcriptional  profiling  experiments,  quantitative  real  timeRT-PCR  analysis  was  carried  out  with  a  subset  of 
the  MMPs  and  ADAMTSthat  characterize  the  expression  patterns  of  each  cluster.  These  data  are  summarized 
Table  4.  As  can  be  seen  in  the  comparison  of  these  data  and  those  presented  in  Figure  5  and  6,  the  expression 
patterns  seen  from  the  prof i  I  i  ng  expert  merits  and  those  determi  ned  by  the  more  accurate  RT-PCR  techni  que 
appear  very  similar.  Theoneexception  was  ADAMTS  1  which  showed  ahigher  level  of  expressi  on  at  day  3  than 
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in  the  profiling  experiment.  These  daa,  in  general,  validate  the  observed  profiles  seen  from  the  clustering  data 
obtained  from  our  transcri  ptiona  profiling  analysis. 

3.  Discussion 

Studies  of  developmental,  pahologica,  and  pharmacological  mechanisms  are  progressively  making  use 
of  transcriptional  profiling  and  computational  analysisto  shed  understanding  on  these  processes  (Su  et  al.,  2002; 
Golub  et  al.,  1999;  Armstrong  et  al.,  2002).  The  biological  complexity  of  fracture  healing  and  its  defined  time 
frame  makes  it  a  particularly  excellent  biological  model  to  apply  transcriptional  profiling.  Two  previous  studies 
have  used  I arge  seal e  transcri pti onal  profi ling  to  study  rat  femur  fracture  heal i ng.  I  n  one  study,  transcri pti onal 
prof i I i was  used  to  identify  how  fracture  healing  might  be  different  in  adult  and  juvenile  rats  The  comparison 
of  the  transcri  pti  onal  profiles  of  the  two  ages  of  rats,  demonstrated  that  aged  rats  had  quantity  vely  lower 
expression  for  a  number  of  bone  marix  mRNAs  including  osteocalcin  and  collagen  type  I.  However,  this  study 
failed  to  identify  changes  in  any  expressed  mRNAs  for  regulatory  factors  that  controlled  the  initial  healing 
process  such  as  BMPs  or  the  early  carti  I  age  formation  during  fracture  healing  was  quantitatively  or  temporally 
different  between  the  two  ages  This  lead  these  authors  to  conclude  that  the  delayed  union  in  olds'  animals  was 
not  due  to  an  altered  repar  process  during  the  period  of  cartilage  formaion,  but  rather  were  related  to  other  genes 
associated  with  the  secondary  bone  formaion  and  remodeling  (Desa  et  al.,  2003). 

A  different  study  by  Hadjiargyrou  et  a.  (2002)  used  suppressive  subtractive  hybridization  cloning  to 
generate  a  cDNA  I  i  brary  of  those  mRNAs  that  are  selectively  up-regulated  i  n  fracture  cal  I  uses  duri  ng  the  ti  me 
course  of  healing.  In  this  study,  cluster  analysis  methods  were  used  to  identify  groupings  of  expressed  mRNAs 
tha  might  be  of  potential  functiona  importance  during  fracture  heaing.  Although  the  full  complement  of  genes 
tha  were  expressed  during  fracture  heaing  was  not  assessed,  severa  nova  gene  families  which  contribute  to 
fracture  heaing  were  identified.  It  is  interesting  tha  using  a  cluster  aralysis  approach  in  this  study,  lead  to  the 
discovery  of  the  involvement  of  the  wnt  family  of  morphogeneti  c  factors  duri  ng  fracture  heal  i  ng.  This 
observaion  isof  particular  interest  given  tha  thisfamily  of  morphogens  has  recently  been  shown  to  be  centra 
reguiaorsof  skaaa  cell  differentiaion  (Boyden  a  a.,  2002;  Bennett  a  a.,  2005;  Glassa  a.,  2005).  These  two 
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studies  then  provide  differing  examples  of  how  transcriptional  profiling  might  be  applied  to  gain  insight  into  the 
biological  mechanisms  that  contributeto  a  complex  processes  such  as  fracture  healing. 

I  n  the  study  reported  here,  a  genome  wide  approach  was  used  to  examine  the  expressions  patterns  of 
mRNAsfor  about  a  third  of  the  known  genes  in  the  mouse  and  a  self  organizing  mapping  approach  was  used  to 
analyze  the  mRNA  expression  patterns.  One  of  the  unique  features  of  our  study  was  that  the  cl uster  analyses 
werevisualizedastwo  dimensional  mosaics  of  colored  tiles  [data  not  shown]  (Echleretal.,2003).  Thistypeof 
analysis  assigns  genesto  tiles  in  a  mosaic  instead  of  the  more  commonly  presented  colored  bars  that  are  used  to 
depict  mRNA  clusters  In  this  w^,  genes  that  are  assigned  to  a  specific  area  in  the  mosaic  based  on  their 
differences  in  expression  relative  to  unfractured  bone,  can  be  easily  be  related  to  those  assigned  in  neighboring 
tiles  Thus,  geneshaving  similar  temporal  and  quantitative  patterns  of  expression  can  be  identified  in  an 
i  ntuiti  vely  obvi  ous  and  easy  to  observe.  Thi  s  approach  gave  us  the  means  of  i  dentifyi  ng  si  multaneous  gene 
activities  acrosstime  without  prior  assumptions  of  any  structure  in  the  data  As  an  example,  wewhereableto 
identify  very  specific  patterns  of  expression  for  several  of  the  MM  PS  and  (MMP2/14/Timp2  and  MMP9/13)  aid 
could  show  that  they  appeared  to  occur  in  atemporally  coordinated  manner.  Concerning  the  expression  profiles 
of  the  ECM  prota  n  mRN  As,  a  di  verse  set  of  mRNAs  encodi  ng  protei  ns  assod  steel  with  basement  membranes, 
connective  tissues  of  blood  vessels,  cartilage  and  type  I  mdrices  were  ail  shown  tofollow  a  similar  temporal 
prof i  I  e  of  i  nducti  on  based  upon  there  positi ons  i  n  the  mosai  c  patterns  of  mRNA  expressi  on.  The  expressi  on  of 
some  of  these  genes  such  as  collagen  types  II  I,  IV,  VI,  XVI,  and  laminin  whileal!  having  been  examined  during 
normal  bone  and  endochondral  development  (Durr  et  al„  1996;  Hausler  et  al„  2002:  Kassner  et  al„  2003;  Luther 
et  al.,  2003),  have  not  been  reported  in  fracture  healing  aid  dearly  now  merit  further  examination. 

Because  the  expression  of  various  M  M  P  mRNAs  showed  such  strong  statistical  correlation  as  a  group  or 
"duster"  that  temporally  moved  together,  these  preliminary  data  directed  the  focus  of  this  study  towards  defining 
the  temporal  relationships  between  the  expression  profiles  of  extracellular  matrix  protei  ns  and  the  larger  group  of 
proteasesand  cellular  fundionsthat  are  involved  in  the  remodeling  and  turnover  of  the  extracellular.  Two 
biological  processesof  fracture  healing  thd  are  critical  to  the  progression  of  the  endochondral  bone  formation  and 
bone  repair  are  tissue  vascularization  and  the  resorption  of  the  endochondral  tissues  (Vu  et  al.,  1998;  Gerber  etal., 
1999;  Zhou  etal.,  2000;  Street  et  al.,  2002;  Ortega  et  al.,  2003).  Previ  ous  studies  have  shown  that  MM  Ps  carry 
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out  essential  roles  during  both  normal  skeletal  development  and  during  adult  fracture  repair  and  are  uniquely 
activated  during  endochondral  development  (Col not  et  al.,  2003;  Holmbeck  et  al.,  2003;  Ortega et  al.,  2003;  Zhou 
et  at.,  2000;  Jmenez  et  al.,  2001).  A  number  of  other  studies  have  also  shown  that  many  of  the  MM Ps  are 
activated  during  the  pathological  turnover  of  cartilage  during  both  rheumatoid  and  osteo-arthritic  processes  (Cho 
et  at.,  2003;  Glasson  et  al.,  2004).  While  many  studies  have  focused  on  individual  MMPs,  only  afew  have  tried  to 
assess  how  the  proteolytic  activities  of  the  various  M  M  Ps  are  coordinated.  Our  studies  would  suggest  that  the 
expressi  on  Tl  M  P2/M  MP2/  MMP14  and  MM  P9/M  M  PI  3  are  two  separate  groups  that  are  coordi  nately  regul  ated. 

I  n  t  hi  s  context  several  recent  studi  es  have  show  n  thdt  M  M  PI  4/M  T- 1 ,  M  M  P2  and  T I M  P2  are  associ  ated  toget  her 
in  a  protein  complex  and  that  pro  MM  P2  is  only  activated  in  conjunction  with  MMP14  and  TIM  P2  (Apte  et  al., 
1997;  Jo  et  al.,  2000;  Karagiannisand  Rope),  2004).  Studies  of  both  pdhological  turnover  of  articular  cartilage, 
as  well  as  normal  developmental  remodeling  of  endochondral  cartilage,  have  shown  that  MMPs  9  and  13  are  al  so 
concurrently  expressed  (Vu  et  al.,  1998;  Lehmann  et  al.,  2005;  Kevorkian  et  al.,  2004). 

Another  finding  from  these  results  is  very  suggestive  to  specific  mechanisms  by  which  calcified  cartilage 
is  resorbed.  These  studies  showed  that  very  low  levels  of  osteoclasts  associated  markers  were  expressed  during 
the  first  ten  days  of  fracture  healing  (Roussel I e  and  Heymann,  2002).  The  observed  elevSion  of  two  separie 
antagonists  to  RANK-L,  osteoprotegerin,  and  TNF  super  family  member  22  was  also  consistent  with  this 
suggesti  on.  TNF  super  f  ami  ly  member  22  i  s  of  consi  derabl  e  i  nterest  and  its  expressi  on  duri  ng  carti  I  age 
remodeling  represents  a  novel  observation,  since  TNFSF22  is  a  decoy  receptor  for  TRAIL  and  can  also  block  the 
activities  of  RANK-L.  The  absence  of  the  expressi  on  of  mRNAsthat  are  associated  with  osteoclast  function, 
concurrent  with  the  increased  expression  of  two  different  antagonist  for  the  key  regulator  of  osteoclast  formation 
(RANK-L)  suggest  that  osteoclast  activity  only  begins  after  initial  periods  of  fracture  hewing  and  cartilage 
f  ormati  on.  Racing  these  results  i  n  temporal  rel  at  i  onshi  p  to  the  expressi  on  data  f  or  M  M  P  supports  the  suggesti  on 
that  the  recruitment  of  osteocl  asts  i  nto  the  carti  I  agi  nous  cal  I  us  ti  ssue  i  s  very  tightly  control  I  ed  and  that  osteod  ast 
activity  only  begins  after  the  initiation  of  the  proteolytic  activities  of  the  MMPs  (Blavier  and  Delaisse,  1995) 

Finally,  the  observations  of  the  reduction  in  the  expression  levels  of  MMP8  also  have  potential  biological 
rel  evanca  SnceMMP8  is  found  predomi  nantly  i  n  neutrophi  I  s,  these  data  woul d  suggest  that  thi  s  cel  I  ty pe  i  s 
largely  absent  in  the  healing  tissues,  until  the  marrow  element  has  been  reestablished  during  the  I  ate  phases  of 
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primary  bone  formation.  MMP8  mRNA  has  also  been  shown  to  be  expressed  at  very  low  levelsduring  the 
pathological  turnover  of  cartilage  associated  with  osteoarthritis  but  is  found  at  higher  levels  in  rheumatoid  arthritis 
suggesting  that  osteoarthritishassomesimilaritiesto  anormal  endochondral  process  (Glasson  et  al„  2004). 
Investigiion  of  the  general  role  of  neutrophils  in  pathological  healing  may  warrant  additional  study. 

Much  less  is  Known  concerning  the  expression  and  mechanisms  of  remodeling  and  resorption  of  cartilage 
proteoglycans  during  endochondral  development,  although  it  is  now  clear  that  members  of  the  ADAMTS  family 
of  proteases  appear  to  be  present  in  cartilage  (Kevorkian  et  at.,  2004;  T ortorella  et  a! .,  1999;  Apte,  2004).  Several 
studies  have  shown  that  that  ADAMTS  4  and  5  are  in  abundance  in  both  in  osteoarthritic  cartilage  and  during 
endochondral  remodeling.  On  the  other  hand,  transgenic  animals  with  deficiencies  in  ADAMTS  4  the  most 
prevalent  of  aggrecanase  in  cartilage,  still  appear  to  undergo  both  normal  development  and  progressive 
osteoarthritic  degradation  under  experimental  conditions  (Glasson  et  al„  2004;  Cal  et  al.,  2002).  The  studies 
reported  here  are  the  first  onesto  datethat  have  examined  the  expression  of  the  ADAMTSgenefamily  during 
fracture  healing  and  showed  like  other  developmental  studies,  that  ADAMTS 4  is  present  during  endochondral 
resorption.  A  second  family  member  ADAMTS  15,  followed  the  expression  of  ADAMTS4.  ADAMTS  15  based 
on  sequence  analysis  would  also  be  considered  an  aggrecanases(Apte,  2004)  although  its  functional  rolein 
endochondral  resorption  is  unproven.  ADAMTS  1  also  has  an  aggrecanase  like  sequences  but  it  reached  maximal 
levelsof  expression  later  than  ADAMTS4  and  15. 

Whi  le  our  initial  experi  ments  were  very  limited  rel  ative  to  the  number  of  ti  me  poi  nts  and  the  f  ami  I  ies  of 
gene  expression  that  were  examined,  this  study  provided  a  convincing  preliminary  set  of  data  related  to  the 
complex  interplay  between  the  ADAMTS,  MMPs  and  re9orptive  activities  of  osteoclasts  during  extracellular 
matrix  remodeling  in  a  fracture  callus  These  initial  data  identified  many  of  the  extracellular  matrix  protein 
mRNAstha  showed  the  largest  quantitative  change  in  their  expression  during  cartilage  and  bone  tissue  formation 
and  suggest  future  lines  of  study  to  examine  specific  interactions  between  cel  Is  types  that  express  various  the 
MMRsor  ADAMTS  and  specific  sdracellular  matrix  proteins  Finally  these  data  demonstrie  how  large-scale 
t ranscri  pti  onal  profi  I  i  ng  data  i  s  abl  e  to  uncover  numerous  compl  ex  rel  ati  onshi  ps  between  vari  ous  mol  ecul  ar 
processes  that  are  occurring  during  fracture  healing.  Future  studies  containing  additional  time  points  and 
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comparing  transgenic  mice  with  defined  deficiencies  in  a  single  gene  will  obtain  more  exact  functional  data 
duri  ng  f  racture  heal  i  ng. 

4.  Experimental  Procedures 

4.1.  Generation  of  murine  tibia  fractures 

Research  was  conducted  in  conformity  with  all  Federal  and  USDA  guidelines,  aswell  as  an  IACUC 
approved  protocol.  All  studieswereperformedonmale8-10weeksoldC57BL/6J(B6)  mice.  Age,  sex  and 
geneti  c  strai  n  of  mouse  were  based  on  previ  ous  data,  i  n  w  hi  ch  the  ti  me  course  of  fracture  heal  i  ng  had  been 
established  by  both  histological  and  selected  candidate  gene  assessment  (Kon  et  al.,  2001;  Cho  et  ai„  2002; 
Gerstenfeld  et  al.,  2003b).  A  modification  of  the  method  of  Bonnarens and  Enhom  (1984),  as  previously 
described  (Kon  et  al.,  2001)  was  used  to  create  si  mpletrans/erse  closed  uni  lateral  fractures  in  the  I  eft  tibia  of  all 
mice.  A  fracture,  whose  configuration  was  inconsistent  with  standardized  placement  criteria  (mid-diaphyseal)  or 
observed  to  be  grossly  comminuted  were  not  used  in  our  studies.  Mice  were  euthanized  at  selected  time  points 
after  fracture  by  CO2  asphyxiation. 

4.2.  Timepointsof  analysis 

Specifi  c  ti  me  poi  nts  over  the  course  of  fracture  heal  i  ng  were  chosen  to  assess  representati  ve  peri  ods  i  n  the 
inflammatory  phase  (day  3);  a  mid  time  point  in  the  endochondral  phase  (da/  10);  and  a  period  of  primary  bone 
formation  and  coupled  remodeling  (day  21).  The  selection  of  these  points  was  based  on  empirical  knowledge  of 
this  model.  The  more  limited  number  of  time  poi  nts  was  chosen  to  amplify  the  andyasto  the  three  peri  ods  when 
specific  biological  processes  would  be  at  their  peak  expression. 

4.3.  RNA  preparation 

After  euthanasia  limbs  were  rapidly  disarticulated  and  surrounding  muscular  was  carefully  removed  with 
care  being  taken  to  not  di  srupt  the  fracture  site  or  cal  I  us  ti  ssuea  Cd  I  us  ti  ssues  were  ci  rcumscri  bed  by  5  mm  on 
either  side  of  the  fracture  aid  an  identical  site  was  excised  from  the  mid  half  of  the  unfractured  contra aterd 
tibiae.  The  fracture  ca  I  uses  were  collected  into  liquid  nitrogen  and  stored  a  -80  °C  until  ready  for  RNA 
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extraction.  Tissueswerepowdered  under  liquid  nitrogen  and  taa  RNA  was  extracted  by  Trizol  (GibcoBRL, 
Gaithersburg,  MD,  US4)  as  previously  described  (Kon  et  al.,  2001).  RNAsto  be  used  for  microarray  were 
redissolved  in  Trizol  and  extracted  a  second  time.  After  isopropanol  precipitation,  RNA  was  dissolved  in  RNase 
freewater  and  further  purified  by  using  RNeasy  mini  columns  (Qiagen  Inc,  Valencia,  CA).  Column  preparations 
were  carried  out  twice  according  to  the  manufacturer’ s  instructions.  At  this  point,  the  RNAs  are  precipitated  with 
Na  Acetate  300mM  and  2.5  volumes  of  100%  EtOH  and  stored  until  use. 

4.4.  Ribonudease  protection  assay  analysis 

For  validati on  of  individual  genes,  ribonudease  protection  assa/s(RPA)  were  carried  out.  All  template 
setswerefrom  PharMingen  Corp.  (San  Diego,  CA).  Assaysfor  bone  extracellular  matrix  genes  were  carried  out 
with  3pg  of  total  mRNA.  Snglestranded  32P- labeled  cRNA  probes  and  denaturing  acrylamide  gels  analysis  were 
as  previously  described  (Choetal.,  2002).  All  data  sets  were  normal  i zed  to  the  housekeeping  gene  L 32. 

4.5.  Reverse  transcriptase  -  PCR 

Identical  RNA  preparations  were  used  for  RT-  PCR  as  for  the  transcriptional  expression  profiling.  All 
reagents  for  the  PCR  analysis  were  from  Applied  Biosystems,  Inc,  and  plate  assays  were  read  on  an  ABI  7700 
Sequence  Detector  (Applied  Biosystems,  Foster  City,  CA).  One  pg  of  total  RNA  was  used  for  each  preparation 
of  cDNA.  All  cDNA  preparations  were  generated  by  random  hexamer  priming.  Reverse  transcription  reactions 
were  initiated  at  25°  C  for  10  minutes  and  subsequently  carried  out  for  an  additional  60  minutes  at  37°  C.  The 
cDNAs  were  diluted  at  a  rati  oof  1 : 50  inRNase-free  water  for  all  PCR  reactions  All  external  and  internal 
pri  mers  w  ere  from  commerci  a  sets  avalablefrom  Applied  Biosystems,  Inc  (Table  1).  Each  plate  contaned  two 
negative  controls  and  a  positive  control  probe.  The  endogenous  control  (p-actin)  wasusedinai  experimental 
samples  for  norma  izai  on  and  quantification  of  each  target  sequence.  All  sampleswere  run  in  triplicates 

Thefractiona  cycle  number  a  which  the  fluorescence  passes  the  fixed  threshold  (CT  vaues)  was  used  for 
quantificaion  by  using  acomparaiveCr  mahod.  This  method  is  described  within  the  Applied  Biosystems,  Inc 
instructiona  manua  for  the  i  nstrument  ai  on.  In  this  mahod,  the  threshold^  vaues)  isdaerminedwhen 
fluorescence  of  probe  reaches  its  exponenti  a  phase  of  accumuiaion.  Sample  vaues  are  then  normaizedtothe 
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threshold  value  for  p-actin  (Actb)  for  each  time-point:  ACT=  X  CT  (exp)  -X  CT  (Actb).  TheCr  value  for  day  0 
was  then  used  asa  reference.  AACT  =  X  CT  (exp)  -  X  CT  (exp  day  0).  The  fold  change  in  mRNA  expression  for 
each  time  point  was  plotted  in  agraph  using  day  0  as  a  reference:  2'AACr(tt3',0)=1. 

4.6.  Microarray  hybridization  and  image  analysis 

Slide  printing,  arra/  labeling,  hybridization,  and  slide  reading  were  performed  at  the  Massachusetts 
General  Hospital  Genomi  cs  Core  Faci  I  ity .  3  i  des  contai  ni  ng  the  Operon  M  us  muscul  us  V 1 . 1  probe  set  (Qi  agen) 
of  over  1 3,000  70-mer  ol  i  gonud  eoti  de  probes  were  pri  rrted  by  the  M  assachusetts  General  Hospital  M  i  croarray 
Core  Fadlity.  All  slides  were  quality  control  tested  after  printing  and  contained  appropriate  positive  and  negative 
control  sequences  for  data  analysis  RNA  from  unfractured  tibiae  was  used  as  the  reference.  Each  experiment 
conta  ned  the  reference  and  a  test  sample.  The  Atlas  Power  script  Fluorescent  Labeling  Kit  (Clontech  Cat  # 
K1860-1)  was  used  to  generate  labeled  cDNAs.  Thetest  and  reference  RNAswere  separately  reverse  transcribed 
into  amine  modified  cDNAs  using  either  Cy3  or  Cy5  labeling  flouros.  9  ides  were  hybridized  simultaneously 
with  both  labeled  probes  and  the  data  was  collected.  A  GenePix  4000B  microarray  scanner  and  software  (Axon 
I  nstruments)  were  used  to  quantify  the  mi  cro  arrays.  Al  I  data  was  stored  as  a  rati  o  of  the  f  I  uorescent  i  ntensity  of 
the  reference  to  the  experi  mental .  The  Bi  oArray  Software  Envi  ronment  (BASE)  (Saal  et  al.,  2002)  was  used  for 
performing  the  quality  control  steps  in  initial  filtering  the  data 

4.7.  Data  normalization  and  filtering 

The  data  was  normalized  by  seal  ingall  individual  intensities  such  that  the  total  intensity  was  the  same  for 
both  comparative  samples  (control  and  fractured)  within  a  single  array  and  across  replicate  arrays  for  that  time- 
point.  The  scaling  is  important  to  eliminate  variation  in  effects  of  differences  in  fluorescent  intensities  with 
different  probes.  Data  from  the  repl  i cate  arrays  were  then  combi  ned  to  i dentify  outl  i ers  and  reduce  statistical 
variation  in  the  data  Genes  with  only  one  data  point  after  initial  selection,  due  to  missing  or  unreadable  data 
were  exd  uded.  The  data  was  combi  ned  usi  ng  the  mean  of  the  repl  i cate  rati  os. 

Using  the  standard  log2  ratio  of  the  two  dye  intensities  from  each  array  spot  representing  each  gene,  the 
distributi  on  of  the  log  rati  os  was  obtained  for  both  the  individual  and  combined  replicates.  Ingeneral  these 
distributions  fall  within  two  SD  within  the  ratio  intensities  plots  and  mean  ratio  and  standard  deviations  are 
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calculated  only  from  these  data  This  allows  one  to  identify  genes  that  were  expressed  sufficiently  above  the  noise 
without  having  to  resort  to  an  arbitrary  mini  mum  ratio  value.  Genes  with  low  intensities  at  the  level  of  the 
background  were  excluded  from  the  analysis.  The  data  was  visualized  using  log  ratio-intensity  versus  log 
intensity  product  on  astandard  R-l  plot,  see  Figure  2  The  plot  revedstheirrtenaty-specificartifadsinthelog2 
ratio  measurements.  The  one  and  two  standard  devi at i on  cutoffs  ai so  reveal  the  i  nherent " noi  se  I  evel "  I i mitati on  i  n 
the  log2  ratio  measurements.  We  carried  out  a  number  of  analyses  for  those  genes  for  which  data  was  available 
across  al  I  ti  me  poi  nts.  Genes  that  di  d  not  change  by  more  than  two  standard  devi  ati  ons  (as  def  i  ned  by  the  Rl  pi  ot) 
in  at  least  onetime  point  were  discarded  to  reduce  complexity  in  later  analysis. 

4.8  Cluster  analvsisof  normalized  data 

Once  a  final  dataset  of  high  quality  probe  measurements  was  attained,  we  used  both  hierarchical  clustering 
(Eisen  et  al.,  1998;  Tamayo  et  al.,  1999)  and  a  self-organizing  maps  (SOM)  (Kohonen  and  Somervuo,  2002) 
based  approach  to  ai  d  i  n  d  usteri  ng  a  mi  I  ar  behavi  ng  genes.  Si  nee  the  results  were  comparabl  e  across  both 
methods,  only  data  from  the  SOM  analysis  is  included  in  this  report.  The  SOM  analysis  was  performed  using  the 
Gene  Express  on  Dynamics  I  nspector(GED  I )  software  package  that  displ^s  high-dimensional  gene  expression 
prof  i  I  es  as  col  ored  coded  mosai  c  i  mages  (Eichl  er  et  al . ,  2003) .  Based  on  positi  on  of  a  genes  pi  acement  withi  n  a 
given  tile  within  the  mosaic,  the  patterns  of  gene  expression  m^  be  easily  identified  without  any  prior 
assumptions  about  the  gene’s  expression.  A  genes  placement  within  a  tile  is  based  on  its  quantitative  level  of 
expression  and  tiles  within  proximity  to  other  tiles  in  the  mosaic  have  similar  quantitative  levels  of  expression. 
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Figure  Legends 

Figure  1.  RNA  expression  anaiya's  of  selected  extracellualr  matrix  protein  mRNAs  across  the  time  course  of 
fracture  healing  by  RPA.  Panel  A)  Representative  autoradiographic  images  of  the  RPA  products  as  resolved  on  a 
6%  PAGE  sequencing  gel.  Days  after  fracture  and  positions  of  each  protected  species  are  denoted.  Panel  B)  The 
graphi  c  analysi  s  of  the  rel  ati  ve  mRNA  I  evel  s  of  sel  ected  genes.  Radi  oacti  ve  counts  of  the  protected  bands  were 
determined  by  direct  p-  counting  on  aflat  surface  beta  detector.  Band  counts  were  normalized  to  the  ratio  of  the 
internal  standard,  L32,  and  expressed  as  a  relative  value. 

Figure  2.  A  representative  example  of  the  spot  analysis  from  one  replicate  for  osteocalcin  expression  is 
presented.  The  ass^r  samples  are  denoted  in  lane  one.  The  raw  fluorescent  intensities thd  were  obtained  with  the 
reciprocal  Cy3/Cy5  cDNA  labeling  and  hybridization  are  shown  in  channels  one  and  two.  The  upper  and  lower 
halves  of  the  read  out  show  the  two  cMa  sets  from  reversing  Cy3/Cy5  hybridizations  between  the  experimental 
and  reference  sampl  es  The  rdi  o  of  reference  channel  to  the  experi  mental  channel  i  s  show  n  i  n  the  thi  rd  I  ane.  The 
descriptive  visual  presentation  of  the  log  variation  form  reference  is  presented  descriptively  as  green  (low)  or  red 
(high).  The  spot  size  and  observed  spots  images  are  shown  in  the  last  two  lanes 

Figure  3.  Rati  o  I  ntensity  pi  ots  of  total  numbers  of  expressed  mRNAs  at  vari  ous  ti  mes  after  fracture.  Rl  pi  ots 
were  generated  as  a  visual  representation  of  uniquely  overexpressed  and  underexpressed  genes  We  used 
unfractured  bone  (day  0)  as  our  reference  and  saw  the  greatest  variance  in  uniquely  expressed  genes  showing  log 
2  fold  greater  or  lesser  from  the  reference  in  the  d^r  3  and  day  10  samples  in  comparison  to  the  day  21  sample. 
The  numbers  of  observed  spots  that  fall  within  two  standard  deviiionsof  conformity  of  normalized  spot  intensity 
are  shown  with  dotted  I  i  nes  1 SD  standard  and  with  a  dashed  I  i  ne  2SD. 

Figure  4.  Graphical  Presentation  of  Various  90Ms  patterns  for  the  MMPs  and  TIMPs.  Four  separate  panels 
showing  the  temporal  patterns  of  expressi  on  of  M  M  P2/M  MP14,  MM  P9/M  M  PI  3,  M  M  P8/M  M  PI  9  and  TI  M  Ps 
1 ,2,and  3  are  shown.  Rdi  o  I  og2  fol  d  change  i  n  expressi  on  rel  ati  ve  to  unfractured  bone  versus  ti  me  after  fracture 
ispresented.  The  groupings  of  individual  MMPSare  denoted  withinthe  keys  contained  in  each  of  the  figures. 
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Figure  5.  Gr^shi  cal  Present  ati  on  of  Various  SOM  s  patterns  for  the  ADAM  Ts.  Three  separate  panel  s  showi  ng  the 
temporal  patterns  of  expression  of  ADAMTS1/10  and  ADAMTS4/15  are  shown  Ratio  log2fold  change  in 
expression  relative  to  unfractured  bone  versus  time  after  fracture  is  presented.  The  groupings  of  individual 
ADAMTS  are  denoted  within  the  keys  contaned  in  each  of  the  figures. 
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Table  1.  RT-PCR  Primers 


Gene 

UniGene  + 

Reference  Number 

Primer  Sequence 

3-Actin 

Mm.297 

NM  007393 

T  ACT  GAGCT  GCGTTTT  ACACCCTTT 

MMP2 

Mm.29564 

NM  008610 

GGACCT  GCAGGGCGGT  GGT  CAT  AGC 

MMP14 

Mm.280175 

NM  008608 

C AGCCCCGA A  GCCT  GGCT GCA  GCA  G 

MMP13 

Mm. 5022 

NM  008607 

AGTT  CCA  AAGGCTACA  ACTT  GTTT  A 

MMP8 

Mm.  1641 5 

NM  008611 

AA  AA  ACT  GCT  GAGA  ATT  ACCTACGA 

MMP9 

Mm. 4406 

NM  013599 

GCAT  CCAGT  ATCT  GT  AT  GGT  CGT  GG 

Adamtsjl 

Mm.  1421 

NM  009621 

GGGA  AAGCCAT  CAGGACCAGGA  AGC 

Adamts 4 

Mm.23156 

NM  172845 

GGACAATGGTTATGGGCACT  GCCT  C 

Adamts  5 

Mm.  112933 

NM  011782 

A  AAAT  ATTACT  CGACAT  CAAGCCAT 

+  The  unigene,  reference  number  and  reverse  primer  sequence  of  genes.  The  forward  primer  and  internal  primer  that  are 
labeled  with  FAM  sequence  are  proprietary. 
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Table  2.  Selective  Expression  Data  From  Micro-Array  Data 

Numeric  data  is  the  Log  2  value  of  expression  relative  to  the  reference  (Un-fractured  bone)* 


Log  2  Over  Expressed 
DAY  3 

3.35  fatty  acid  binding  protein  5,  Fabp5 
3.02  follistdtin-likeFstl 

272  osteoblast  specific  factor  2  (faeddin  l-like)  OSF-2 

2.45  TNF  receptor  superfamily  member  22 

0.52  TNF  receptor  superfamily,  member  11b  (OPG) 

DAY  10 

4.17  cartilage  oligomeric  matrix  protein  (Comp) 

4.05  osteoglydnOgn 

3.82  mdrix  gamma-carboxygl  utamate  (gl  a)  protei  n 
3.7  fol  I  istati  n-l  i  ke  Fstl  hepari  n  bi  ndi  ng 

2.98  pelecan  (heparan  sulfde  proteoglycan  2)  Flspg2 
2.79  midkine  Mdk  growth  factor;  heparin  binding 
2.65  insulin-li  ke  growth  factor  bi  ndi  ng  protei  n  5 1  gf  bp5 

1.63  TNF  receptor  family  superfamily  member  22 
0.79  TNF  receptor  superfamily,  member  11b(OF>G) 

DAY  21 

2.26  twist  gene  homolog  Twist  DNA  binding;  protein 
2.25  osteoglydn  Ogn 

1.02  TNF  receptor  family  superfamily  member  22 
0.52  TNF  receptor  superfamily,  member  llb(OFG) 


Log  2  Under  Expressed 


-2.16  bone  gamma-carboxygl  utamdte  protein,  Osteocalan 
-218  acid  phosphdase  5,  tartrate  resistant  Acp5  hydrolase 
-3.30  TATA  box  binding  protein 
-4.44  carbonicanhydrase2Car2lyase 


-270  FGF  indudble  16  Fin  16 
-2.72  high  mobility  group  box  2  Flmgb2 
-5.69  carboni  c  anhydrase  2  Car2 1  yase 


-262  neutrophi  lie  granule  protein  Ngp 
-3.18  eosinophil-assodated  ribonudease3  Ear3 


*  All  values  represent  theincfividual  mRNAs  express)  on  as  a  ratio  between  unfractured  bone  (reference)  and  the 
expei mental  sample  All  valuesarealso  presented  asthelog  2  of  thevalue  As  an  example  COM  Pis  expressed  as  4. 17 
log2  fold  greyer  at  day  ten  or  ~21  fold  higher  in  expression  than  unfractured  bone  If  the  mRNA  is  unde  expressed  the 
numeica!  value  is  log2  fold  I  ess  than  urrfradured  bone  and  i  s  e<  pressed  as  a  negdi  ve  vai  ue. 
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Table  3.  Prevalent  GEDI  Groupings  Showing  Elevated  Expression  Throughout  The  Repair  Period 


A.  Group  I 

Generalized  ECM  Proteins 

Procollagen  typeV  alpha  1 

Procollagen  type  V  alpha  3 

Procollagen  type  VI  alpha  2 

Procollagen  type  VI  alpha  1 

Procollagen  type  VI,  alpha  3 

Procollagen  type  XI ,  alpha  1 

Procollagen  type  XVI,  alpha  11 

Procollagen-C  endopeptidase  enhancer 

Collagen  triple  helix  containing  1 

Procollagen-Lysi  ne  2-Oxogl  utarate  5-Di  oxygenase 

Osteoglydn 

Cartilage  Associated 

Cartilage  oligomeric  matrix  protein 

Perlecan 

procollagen  gene  (IX)alpha2 
procollagen  type  (IX)alphal 
Ecml :  Extracellular  matrix  protein  1 

Non  Defined  Families 
Osteoblast  specific  factor  2 
S100  Calcium  binding  protein 
Cadherin  4 
Vimentin 
Annsdn  A6 


B.  Group  II 

Procol  lagen  type  1  al  pha  1 
Procollagen  type  I,  alpha2 
Procollagen,  typeV,  alpha 2 
Procollagen  C-proteinase  enhancer  protein 

Vessel  and  Basement  Membrane  Related 

Procollagen  type  IV  alpha  1 

Procollagen  type  IV,  alpha2 

Thrombospondin  4 

Thrombosjxmdin  2 

Thrombospondin  1 

Laminin  beta 2 

Lam  nin  gamma  1 

Lamin  A 

Nidogen  1 

Procol  I  agen  type  1 1 1  al  phal 
Fibrillin 

TNF  Family 

TNF  receptor  superfamily  member  22  (TRAIL  Decoy) 
TNF  receptor  superfamily,  member  11b 
(osteoprotegerin) 
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Table  4.  RT-PCR  Fold  Induction  Over  Time 


A  Time 

Gene 

MMP2 

MMP8 

MMP9 

MMP13 

MMP14 

Adamtsl 

Adamts4 

Adamts5 

DAYO 

1 

1 

1 

1 

1 

1 

1 

1 

DAY  3 

8.397  + 

0.1 

0.2 

0.3 

6.06 

5.76 

15.24 

3.6 

DAY  10 

19.159 

0.3 

4.25 

10.1 

10.126 

7.889 

19.42 

3.45 

DAY  21 

5.098 

0.645 

4.8 

7.315 

4.99 

5.57 

7.67 

3.317 

1  The  fold  change  in  mRNA  as  normalized  to  fl-actin  and  expressed  as  a  fold  change  relative  to  unfractured  hone 
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MICROABSTRACT:  The  spatial  temporal  cellular  activities  of  fracture  callus  morphogenesis  were 
defined  using  serial  sections,  standard  histological  measurements,  in  situ  hybridization  and  three 
dimensional  tissue  reconstruction.  These  data  showed  that  endochondral  bone  formation  occurred  in  an 
asymmetric  manner  within  femur  fracture  calluses  such  that  cartilage  tissues  formed  predominantly 
distal  and  lateral  to  the  fracture  and  that  the  spatial  confines  of  the  cartilage  formation  established  the 
spatial  patterns  of  each  of  the  subsequent  stages  of  bone  healing. 

Introduction:  Tissue  morphogenesis  is  control  led  by  spatial  temporal  relationships  between  various  cel  Is  types 
and  their  secreted  morphogena  Within  a  growth  piatethe  spatial  arrangements  between  the  various  cell  types 
that  contri  bute  to  the  skel  etal  morphogenesi  s  are  defi  ned  i  n  I  arge  part  by  the  ti  ssue’ s  I  i  near  pattern  of  growth. 

Whi  I  e  fracture  repai  r  recapitul  ates  many  aspects  of  endochondral  bone  formati  on  the  morphogeneti  c  pattern  of 
the  various  tissues  that  contri  bute  to  a  fracture  healing  are  not  as  easily  discernable. 

Methods:  Serial  histological  sections  of  mouse  and  rat  fracture  cal  I  uses  were  used  to  define  the  spatial  and 
temporal  changes  in  tissue  composition  associated  with  the  progression  of  fracture  healing.  In  situ  hybridization 
to  collal  and  col2a1  mRNAs  were  used  to  define  areas  of  new  cartilage  and  bone  formati  on.  Two  dimensional 
representations  of  tissue  composition  were  compared  to  three  dimensional  computer  assisted  reconstructions 
using  the  same  serial  sections  Reconstructed  volumes  were  used  to  define  the  spatial  development  of  the 
morphogeneti  c  f  i  el  ds  thd  gave  ri  se  to  the  fracture  cal  I  us. 

Results:  Both  compositional  and  three  dimensional  reconstructions  demonstrated  that  cartilage  tissue 
developed,  in  an  asymmetrical  pattern  around  the  fracture.  After  14  days  of  fracture  healing,  approximately 
70%  of  the  volume  of  callusthat  formed  in  a  mid-diaphyseal  femoral  fracture  was  on  the  distal,  lateral  aspect  of 
the  bone.  As  the  cal  I  us  devel  oped  proxi  mal  ly ,  the  carti  I  age  tissue  rotted  toward  the  medi  al  bone  surface. 
Subsequently  as  fracture  healing  progressed,  trabecular  bone  replaced  the  calcified  cartilage  proceeding  from  the 
edges  of  the  cal  I  us  i  nward  towards  the  fracture  with  the  area  around  the  fracture  to  be  the  I  ast  to  be  repl  aced. 

The  asymmetric  pattern  of  endochondral  bone  formation  was  confirmed  by  reconstructing  the  three  dimensional 
pattern  of  collagen  type  I  (collal)  and  II  (col2A1)  mRNA  expression,  as  determined  by  in  situ  hybridization  in 
mouse  femur  fractures.  These  results  corroborated  those  in  the  rat  showing  the  same  a^mmetric  pattern  of 
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endochondral  bone  f  orm^i  on.  Three  di  mensi  onal  reconstruct  on  of  the  cal  I  uses  subsequent  to  the  resorpti  on  of 
the  cartilage  showed  a  complex  structure  in  which  thin  shell  of  new  exterior  cortical  bone  was  now  supported  by 
an  i  nteri  or  structure  of  trabecul  ar  bone  that  had  repl  aced  t  he  ori  gi  nal  carti  I  age  ti  ssues. 

Conclusions:  These  results  demonstrate  that  specific  developmental  mechanisms  control  the  asymmetric 
pattern  of  cal  I  us  formation  and  suggest  that  this  pattern  recapitulates  the  original  asymmetry  of  endochondral 
tissue  development  of  the  femur.  On  the  other  hand  the  spatially  coupled  mechanisms  that  remodel  the  callus 
from  a  greater  cross  sectional  diameter  to  a  lesser  diameter  must  be  different  from  the  developmental 
mechanisms  that  are  observed  during  appositional  bone  growth. 

KEY  WORDS: 

Orthopaedics,  Fracture  Repair,  BoneHistomorphometry,  Growth  and  Development 
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INTRODUCTION 


Bone  i  s  one  of  the  few  ti  ssues  i  n  the  body  that  undergoes  true  regenersti  on  i  n  response  to  injury  and 
many  of  the  mechanisms  involved  in  skeletal  repair  appear  to  recapitul  ate  the  events  of  embryo!  ogic 
developmental.12^  Our  present  understanding  of  the  basic  mechanisms  of  skeletal  repair  suggests  that  itisa 
multifaceted  process  involving  complex  interactions  between  the  immune,  hematopoietic,  and  vascul a 
systems4'5'6'7'8  These  systems  all  interact  with  the  underlying  mesenchymal  cell  population  within  and 
surrounding  bone  and  contribute  to  theinduction,  growth,  and/or  maintenance  of  the  skeletagenic  cellular 
responsesthat  regenerate  the  injured  tissue.  The  overall  spatial  relationships  between  the  morphogenetic  fields 
and  the  various  proximate  tissuesthat  contribute  to  the  repair  process  (Mine  the  final  tissue  geometry  aid 
composition  of  the  fracture  callus  Theend  result  of  these  developmental  processes  is  the  reestablishment  of 
both  the  original  geometry  and  biomechanical  competency  of  the  skeleton. 

Because  the  size,  shape  and  mataial  properties  of  adult  bones  are  determined  ealy  in  life  and  ae  in  pat 
controlled  by  intrinsic  genetic  factors, 9'iai1  it  isimportant  to  understand  the  biological  processes  (or  factors)  that 
contri  bute  to  the  vari ati  on  i  n  these  morphometri c  and  maeri  a  traits  Thi  s  can  be  accompl  i  shed  i  n  the  context  of 
the  devaopmenta  processestha  ae  acti vaed  during  fracture  repar.  A  quantitaive definition  of  the  spaia 
geomary  and  naure  of  the  tissue  formai  on  (endochondra  versus  osseous)  as  well  as  a  definition  of  the  spaia 
naure  of  tissue  remodeling  tha  regenaaesthe  origina  bone  is  a  so  necessary  for  our  understanding  of  fracture 
heaing.  In  this  study,  we  conducted  a  histomorphometric  anaysisof  thetempora  progression  of  fracture 
heaing.  We  defined  tissue  compositiona  measurements  the  statical  accuracy  of  these  measurements  aid 
made  comparisons  of  these  measurements  to  three  dimensiona  histologica  reconstruction  of  fracture  ca  I  us 
ti  ssues  a  vari  ous  devel  opmenta  stages  duri  ng  fracture  heal  i  ng. 

In  order  to  link  histologic  changes  with  the  expression  of  marix  genes  critical  to  the  repar  process  the 
tempora  and  spaia  expression  of  collagen  types  I  and  II  mRNAs  were  assessed  by  in  situ  hybridizaion.  The 
results  provide  the  basic  anaytica  bridge  between  the  cellular  and  molecular  processestha  form  and  mantan 
skaaa  tissues  and  the  geometric  and  maeria  naure  of  theca  I  us  tissue  tha  produces  its  mechanica  properties. 
The  <Ma  further  demonstraetha  the  spaia  arrangement  of  varioustissues  within  the  cal  I  us  develop  aid  hea  in 
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an  asymmetric  pattern.  This  asymmetry  is  suggestive  of  the  original  spatial  pattern  of  growth  in  the  uninjured 
long  boneandm^  correlate  with  the  activity  of  the  proximal  or  distal  physes. 


MATERIALS  AND  METHODS 

Production  of  Simple  Transverse  Fractures 

Animal  research  was  conducted  in  conformity  with  all  federal  and  USDA  guidelines,  as  well  as  an 
IACUC  approved  protocol.  MaleSprague-Dawley  rats(n  =  3  per  group)  we'ghing  449±39g  and  about  7-9 
months  old  were  used.  In  experiments  assessing  mice,  mal  e  C57BL/6J  mi  ce  (n  =  5  per  group)  at  8to  lOweeks 
of  age  were  used.  All  ratswere  purchased  from  Harlan  Bioproducts  for  Science,  Indianapolis,  IN  and  all  mice 
were  purchased  from  Jackson  Laboratories,  Bar  Harbor  ME  Closed,  simple,  mid-diaphyseal,  transverse 
fractures  of  the  rat  femurs  were  produced  as  described  by  Bonnarens  and  Einhorn.12  A  modification  of  this 
procedure  was  used  for  the  mice5.  However  unlike  our  previous  studies  in  mice  in  which  the  tibia  was 
fractured,  the  femur  was  fractured  i  n  the  current  studi  es.  I  n  the  case  of  the  mouse  femur  fractures, 
intramedullary  pi  ns  were  inserted  retrograde  through  the  distal  condyle  of  the  femur  and  stabilized  in  the  major 
trochanter.  Radiographic  assessment  of  the  fractures  was  performed  immediately  after  fracture  and  at  the  time 
of  euthanasia  Fractureswhich  did  not  occur  in  the  mid-diaphysis  or  were  excessively  comminuted  were  not 
used  in  the  study.  All  ani mal s  were  euthani zed  by  CO2  asphyxiation  at  7, 14, 21  and  35  d^rs  post-fracture.  At 
the  ti  me  of  euthanasi  a  fracture  cal  I  us  sped  mens  were  harvested  and  carefully  d  eaned  of  musd  e  and  soft 
connedivetissue.  Cal  I  us  dimensions  were  measured  at  this  time  in  anterior- posterior  (AP)  and  medial -lateral 
(ML)  dimensions  using  an  electronic  caliper  (#R163-7251006C,  Maryland  and  Metrics*  Baltimore,  MD). 


Tissue  Fixation  and  Decalcification 

For  histological  assessments,  the  bones  with  a  small  amount  of  surrounding  musd  e  and  soft  tissues  were 
fixed  for  2  days  in  4%  paraformaldehyde  in  phosphate  buffered  saline  st  4°C.  Sped  mens  were  completely 
decaldfied  in  14%  w/v  EDTA  changed  three  times  per  week  for  approximdely  2-3  weeks  while  shaking  at  4°C. 
$>ed  mens  were  rinsed  and  placed  in  ice  cold  PBS  for  final  processing.  After  decal  dfi  cation,  theintranedullary 
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pin  was  removed  and  the  anatomic  center  of  the  callus  was  determined  from  the  x-ray  measurements  and  was  set 
as  the  poi  nt  of  the  transverse  fracture.  The  rat  bones  were  cut  transversely  with  a  sharp  seal  pel  at  two  poi  nts 
5mm  proximal  and  5mm  distal  to  the  center  of  the  fracture  callus  Another  cut  was  then  made  at  the  center 
where  the  original  fracture  had  been  determined  from  the  X-ray  assessments  thereby  creating  two  half  callus 
sped  mens  M  ouse  cal  I  uses  were  i  sol  ated  i  n  the  same  manner  but  were  I  eft  i  ntact.  After  dehy  drati  on  i  n  graded 
concentrations  of  ethanol  to  100%,  the  sped  mens  were  transferred  to  xylene,  and  embedded  in  a  step-wise 
manner  under  vacuum  in  50%  xylene  paraffin,  then  100%  paraffin.  The  two  halves  of  the  rat  callus  were 
positioned  in  a  angle  block  of  low-melt  paraffin  with  the  fracture  centers  fad ng  the  cutting  surface.  The  mouse 
cal  I  uses  were  embedded  intact  with  the  distal  end  (end  dosest  to  thetibial  joint)  fadng  the  cutting  surface  A 
counting  microtome  was  used  such  that  the  total  linear  distance  thd  had  been  sectioned  through  a  block  could  be 
monitored.  Five  micron  thick  paraffin  sections  were  cut  and  placed  on  poly  L-lysine  coated  slides,  dried 
overnight  and  used  immediately  or  stored  at  4°C. 

Histomorphometric  Parameters  and  Sampling  Scheme 

Three  rat  or  five  mouse  bones  were  used  for  each  time  point  of  the  histomorphometric  analysis.  Twelve 
to  fifteen,  5  ^m-thick  serial  sections  were  taken  at  100pm  increments  across  the  entire  length  of  each  callus 
(~6000^im).  In  this  manner  a  series  of  -60  segments  containing  groups  of  serial  sectionswere  obtained  that 
spanned  the  complde  length  of  the  callus  The  positions  of  each  segment  from  which  serial  sections  wastdren 
is  denoted  relative  to  the  center  (C)  which  is  the  site  of  the  fracture  and  various  segments  containing  groups  of 
serial  sections  are  enumerated  in  distances  proximal  (F)  or  distal  (D)  to  the  center  as  defined  by  the  original 
anatomy  of  the  femur.  The  hi  stomorphometri  c  sampl  i  ng  schemes  that  were  used  to  exami  ne  cal  I  us 
morphogenesis  and  the  progression  of  healing  are  summarized  in  Figure  1 A  and  B.  The  histomorphometric 
indicesthat  were  examined  were  as  previously  described13  and  included:  callus  diameter  (CDm),  callus  areas 
(C/Ar)  cartilage  areas  (Cg/Ar),  void  areas  (V/Ar)  [inclusive of  empty  space  and  hematopoietic  elements],  total 
osseous  tissueTOT/Ar  [inclusive  of  the  original  cortical  bone  thd  was  present  prior  to  the  fracture  new  cortical 
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bone,  trabecular  bone  and  lining  cells],  and  osteoclast  density  (Oc/Ar).  The  abbreviations,  nomenclature  and 
unit  measurements  are  in  accordance  with  the  currently  used  standardsfor  histomorphometry  of  intact  bone.14,15 

Staining  Methods  and  Image  Collection 

Sections  were  stained  with  Safranin  O-Fast  Green  as  previously  described16  to  discriminate  mdure 
cartilage  from  bone  and  non-cartilage  connective  tissues.17  Tartrate  resistant  add  phosphatase  (TRAF),  a  marker 
for  osteoclasts,  was  detected  using  an  azo-dye  coupling  method  with  the  slight  modification  of  fast  red  violet  LB 
salt  (Sigma  F-3381,  Sigma  Chemicals,  St.  Louis,  MO)  repladngthefast  red  TR  salt.  Each  section  was 
photographed  with  an  Olympus  BX51  light  microscope  attached  to  a  digital  camera  at  10.25X  and  downloaded 
into  an  I  mage-Pro  Plus  Version  4.1 .0.0  for  Windows  program.  An  area  of  interest  (AOI)  was  created  by 
loading  a  uniform  box  (5mm  x  7.9mm)  onto  the  photograph  and  centering  thecallus  within  the  box.  The  bone 
wasthen  outlined  within  the  uniform  AOI,  excluding  any  muscle,  soft  tissue  or  periosteum.  Usingacolor 
match  program,  the  total  area  of  cartilage  (red)  and  bone  (green)  was  initially  identified  quantified  using  afilter 
range  of  573.9p.m-5.74e<-03jim.  Subsequently  colored  matched  areas  defining  specific  tissues  were  individually 
assessed  and  final  areas  were  hand  traced.  Mean  values  were  calculated  for  sped  mens  sampled  from  each  bone 
and  then  used  to  create  group  means,  standard  deviations  and  standard  errors  both  with  resped  to  time  after 
fracture  and  animal  group. 

For  measuring  either  TRAP-stained  cells,  contiguous  areas  were  sampled  colleding  as  many  as  7 
micrographs  for  each  transverse  sedion,  such  that  when  overlapped  with  each  other  they  spanned  the  width  of 
thecallus  (Figure  1C).  Microscopic  images  were  captured  as  described  above  and  adjacent  microgr^hs 
(10.5mm  x  7.8mm)  were  taken  using  the20X  objedive  and  downloaded  onto  an  Image-Pro  PI  us  Version  4. 1.0.0 
for  Windows  program.  Using  acolor  match  program,  the  osteodasts  from  each  photograph  were  quantified 
ua  ng  a  spot  fi Iter  and  a  pre-defi  ned  pixel  val  ue  that  def i  ned  the  lower  and  upper  di mensi ons for  the  spots 
counted.  Mean  values  of  TRAP  positive  cells  per  callus  age  groups  were  calculated  from  cross  sedional 
measurements  taken  at  1000  micron  increments  across  the  total  longitudinal  distance  of  each  cdlus(N=37-49 
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mi  crographs  per  cal  I  us) .  The  numbers  f  rom  the  multiple  i  mages  per  group  were  combi  rted  and  sped  men  means 
calculated  and  used  to  crede  group  means  standard  deviations  and  standard  errors 

Statistical  Methods 

Each  callus  was  divided  into  approximately  30  segments  proximal  and  30  segments  distal  with  each 
segment  of  serial  sections  taken  at  1 00  pm  intervals  from  the  anatomic  center  as  defined  by  the  fracture.  Thisis 
shown  schematically  in  Figure  1A.  Each  segment  contained  a  group  of  12-15  serial  5pm  sections.  One  section 
from  each  segment  was  then  analyzed  such  that  the  entire  length  of  each  callus  (total  of  -60  sections  one  from 
each  segment  acrossthe  entire  length  of  the  callus)  were  used  for  quantitative  measurements.  Thestatisticd 
vari  abi  I  i  ty  i  n  the  van  ous  hi  stol  ogi  cal  measurements  i  n  each  cal  I  us  was  obtai  ned  by  determi  ni  ng  the  meai  val  ues 
from  the  measurements  made  from  the  -  60  sections  taken  acrossthe  entire  length  of  each  of  the  calluses  This 
was  then  used  to  defi  net  he  so- cal  l  ed  reference  val  ue  for  each  of  the  measurements  The  95%  percent 
corif  i  dence  i  rrterval  was  calculated  around  each  reference  val  ue  and  was  used  to  defi  ne  the  mi  ni  mat  number  of 
sections  per  cal  I  us  to  accurately  predict  val  ues  for  each  measurement.  The  minimal  number  of  segments  from 
which  sections  should  betaken  was  estimated  by  taking  the  sections  from  increasing  increments  (every  2nd,  3rd 
4th  etc  segments)  starting  with  the  central  proximal  (PI)  and  distal  (D1)  sections  This  is  presented 
diagrammatically  in  figure  1  A,  with  variations  in  the  sampling  measurements  denoted  at  progressively  larger 
i  ncrements  withi  n  the  segments  denoted  by  X&  The  mean  val  ues  for  each  of  the  pri  may  aiatomi  cal 
measurements  were  then  computed  for  each  set  of  sample  and  compaed  to  the  95%  confi dence  i  ntervd  for  the 
respective  reference  value.  From  these  preliminary  statisti cal  analyses  it  was  determined  that  sampling  st  ai 
increment  of  not  greyer  than  500  pm  and  spaced  uniformly  across  the  complete  longitudinal  distance  of  the 
cal  I  us  accurately  fell  withi  n  the  95%  confi  dence  i  nterval  of  the  references  for  each  of  our  measurements.  The 
means  and  standad  deviations  between  all  histological  measurements  in  the  daa  presented  in  this  study 
thaefore  were  all  derived  from  a  sampling  increment  of  500  pm  uniformly  spaced  acrossthe  length  of  the 
callusi 
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It  should  be  noted  as  well  that  other  sampling  schemes  such  as  taking  sectionsjust  from  segments 
around  the  center  for  example  do  not  provide  means  that  fit  with  the  confidence  interval  with  as  few  a  number  of 
sections  as  uniformly  spacing  the  increment  across  the  entire  length  of  the  cal  I  us.  Fewer  analysis  for  carrying 
out  futwe  studies  were  also  determined  (data  not  shown). 

In  Situ  Hybridization 

Five  micron-thick  sectionswere  placed  on  poly  L-lysine  or  15%  Bond-fast  glue  (Manco  Inc.,  Avon, 

OH)  coded  slides,  dried  overnight  and  used  immediately  or  stored  at  4  C.  Both  antisense  and  sense  control,  ^S- 
labeled,  cRNA  probes  were  used  for  hybridization.  Linear  cDNA  sequences  were  incubated  with  either  T7  or 
SF6  RNA  polymerase  in  the  presence  of  35S-UTP(New  England  Nuclear,  Boston,  MA),  united  nucleotides, 
10mM  DTT,  and  Rnasin  RNase  inhibitor  (Promega,  Madison,  Wl).  Labeled  cRNA  probes  were  purified  using 
Mini  Quick  Spin  RNA  columns  (Roche  Diagnostics,  Indianapolis,  CA).  Prehybridization:  9  ides  were 
deparaffinized  in  xylenes  followed  by  rehydration  in  graded  ethanol  solutions,  rinsed  in  0.85%  NaCI  (5 
minutes)  and  IX  PBS  (5  min),  and  treated  with  proteinase  K  (20pg/ml)  for  8  minutes  at  37°C.  SI  ides  were 
dipped  successively  in  IX  PBS  (5  min),  4%  paraformaldehyde  (5  minutes),  acetylated  in  0.25%  acetic  anhydride 
in  100  mM  triethanolamine  (10  minutes),  washed  in  IX  PBS (5  minutes),  0.85%  NaCI  (5  minutes),  dehydrated 
in  graded  dhanol  solutions,  and  air  dried.  Hybridization:  Hybridizdion  solution  contains  50%  deionized 
formamide,  0.3M  NaCL,  20mM  Tris-HCI  (pH  7.5),  5mM  Na^EDTA  (pH  8.0),  10%  dextran  sulfate,  IX 
Denhardt’ s  solution,  0.5  mg/ml  total  yeast  RNA,  and  lOOmM  DTT.  ^labeled  cRNA  probes  Ere heated  at  80°C 
for  two  minutes  and  placed  on  ice.  Hybridization  solution  containing  ^labeled  probe  at  Sx^cpm/pl  was 
added  to  each  slide  in  a  50  -  60pl  volume,  a  plastic  HybriSip  (Research  Products  International,  Mt.  Prospect, 

I L,  USA)  was  placed  on  each  slide,  and  the  si  ides  were  incubated  at  52°C  for  16  hours  After  hybridization, 
si  ides  were  washed  in  5X  SSC  for  30  min  at  50°C,  2X  SBC  with  50%  deionized  formanidefor  20  min  st  65°C, 
and  rinsed  twice  in  10X  STE,  lOmin  each,  at  37°C.  Sides  were  then  treated  with  a20pg/ml  solution  of  RNase  A 
(Pharmacia  Biotech,  R  scat  aw  ay,  N^>  for  30  minutes  at  37°C,  washed  in  2X  and  0.1X  SSC  for  five  minutes  each, 
dehydrated  in  graded  ethanol,  and  air  dried.  Autoradiography:  Air  dried  si  ides  were  dipped  in  Kodak  NTB-2 
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emulsion  (Eastman  Kodak,  Rochester,  NY),  drained,  air-dried  for  1  hour,  and  placed  in  alight-proof  container 
with  desiccant  a  4°C  for  2-3  weeks.  3  ides  were  da/eloped  in  Kodak  D-19  developer,  fixed  in  Kodak  fixer,  and 
count erstained  with  Fast  Green  and  Safranin-O. 

Three  Dimensional  Reconstructions 

Three  di  mensional  reconstructions  were  performed  usi  ng  the  Ami  ra  3.0  software  system  for  3D 
visualization  data  analysis  and  3D  geometric  reconstruction  (Visual  Concepts  GmbH,  Company,  Konrad-Zuse- 
Zentrum  (ZIB),  Research  Institute  Berlin,  Germany).  All  image  modifications  were  performed  in  Adobe 
Photoshop  Software,  (Adobe  Systems,  Incorporated;  San  Jose,  CA  USA).  In  order  to  use  the  col  or  recognition 
operations  employed  during  the  3D  reconstruction  process,  soft  tissue  a’eas  surrounding  the  exteriors  of  the 
calluses  were  digitally  removed  from  incfi  vidud  images  by  hand  tracing  tissue  perimeters  in  the  Photoshop 
program,  and  cutting  these  areas  from  the  images  Amira  3.0  3D  reconstruction  software  is  unableto  recognize 
multi-channel  images  that  are  acquired  initially,  but  subsequent  grayscale  outputs  yield  enough  contrast  to 
differentiate  the  tissues  of  interest.  In  order  to  maximize  the  color  contrast  differentials  within  a  gray  scale 
image,  the  color  differentials  were  first  digitally  extracted  to  optimize  staining  contrasts  Since  areas  of  mdure 
cartilage  tissue  stained  bright  red  with  SafraninO/  Fast  Green  areas  were  hand  traced  and  the  green  channel  was 
fi  rst  extracted  from  the  pri  mary  RGB  i  mage  thereby  enabl  ing  the  carti  I  age  to  be  easily  contrasted  as  black  areas 
i  n  separate  gray  seal  e  i  mages  Thi  s  hi  gh  degree  of  contrast  i  n  the  green  channel  makes  this  staining  property 
ideal  for  accurate  measurements  using  imagethresholcfing  techniques.  For  cortical  bone  measurements  aid 
image  discriminations  these  aeas  were  outlined  in  each  si  ice  using  a  manual  segmentation  tool.  Once 
individual  imageswere  modified  they  were  imported  into  Amira  software  program  and  the  slices  manually 
al  i  gned  usi  ng  edge  to  edge  matchi  ng  of  the  edges  from  the  surgi  cal  pi  n  hoi  e  i  n  the  marrew  cavity  using  Amira’s 
si  ice  al  i  gnment  editor.  The  si  i  ce  a  i  gner  permits  user  vi  sual  izati  on  of  two  consecutive  si  i  ces  at  a  ti  me,  each 
exhibiting  opposite  contrast.  Thisyieldsa50%  grey  image  representing  two  perfectly  mdched,  overtyping 
images.  The  editor  a  lows  for  berth  dices  to  be  aigned  with  ether  the  previous  or  the  next  si  ice  as  well  as  the 
ability  to  mani  pul  ae  the  entire  stack  of  dices  in  order  to  center  the  field  of  view  within  a  bounding  box. 
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In  order  to  visualize  3D  surf  aces,  the  different  tissues  were  labeled  or ‘segmented’  using  the 
segmentation  editor  in  the  software  program.  Individual  sets  of  I  dads  were  created  for  cartilage,  cortical  bone, 
the  cal  I  us  as  a  whol e  and  for  the  space  between  the  trabecul ae.  The  segmentati on  editor  allows  for  both  manud 
and  automated  segmentation.  For  most  tissues  it  was  possible  to  use  the  automated  segmentation  to  label  slices 
usi  ng  threshoi  di  ng.  Threshol  di  ng  provi  des  a  si  mpl  e  and  fast  method  of  I  abel  i  ng  without  fal  se  positi  ve  results 
False  positives  were  easily  identified  by  comparing  the  labeled  images  with  the  originals  Thenon-resampled 
stack  of  si  i  ces  i  s  used  for  this,  and  the  labels  are  I  der  resampl  ed  to  maxi  mi  ze  accura^ . 

The3D  spatial  distribution  of  collagen  type  I  -  and  col lagen  type  II  -  expressing  cells  in  the  healing 
fracture  calluseswere  visualized  using  in  situ  hybrid  zati  on.  Tosimplify  the  reconstructions  of  thecoll^en  I 
and  collagen  1 1  in  situ  hybridizations,  a  color  select  was  performed  for  these  markers  in  Adobe  Photoshop  5.5 
usi  ng  dark  f  i  el  d  i  mages  As  with  the  col  or  stai  ni  ng,  thi  s  pre-opt  i  mi  zati  on  of  the  i  mages  was  used  to  enhance  the 
color  differentials  in  the  gra/scaie  i  mages.  Through  thresholding  a  binary  image  is  produced,  thereby 
simplifying  the  segmentati  on  process  The  separate  stacks  of  si  ices  for  col  lagen  I  and  collagen  II  in  situ 
hybri  di  zati  on  i  mages  were  i  nserted  i  nto  the  exi  stent  stacks  generated  from  the  three  render!  ng  of  the  I  i  ght 
images,  paring  careful  attention  to  the  spatial  coordinates  to  ensure  a  proper  overlap. 

Volumes  were  obtained  from  each  reconstructed  surface  using  an  automated  volume  measurement 
function  in  Amira  The  scale  was  determined  by  comparing  the  known  width  of  each  im^e  with  the 
measurement  determined  by  Arrtrafor  the  width  of  the  bounding  box.  A  stage  micrometer  was  used  to 
determine  the  actual  width  of  each  image  in  Image-Pro  Plus  The  scale  was  thus  determined  as  being  100pm  = 

1 .32  Amira  Units,  and  hence  1  Amira  Unit3  =  0.000435  mm3 


RESULTS 

Quantitative  Analysis  of  the  Spatial  Morphogenesis  of  Callus  Development 

The  cross  sectional  diameters  and  area  measurements  of  the  calluses  reflect  the  remodeling  of  the 
tissues  asthty  diminish  from  their  peak  dimensions  towards  their  original  anatomical  sizes  and  configurations 
(Figure2).  Fortherat,  maximal  fracture  cal  I  us  tissue  areas  were  deserved  at  21  c%s  after  which  their  size 
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progressively  diminshed  over  the  time-course  of  healing.  It  is  interesting  to  notethat  while  the  medid-latera 
di  ameters  showed  the  greatest  changes  over  the  heal  i  ng  peri  od,  the  posterior-anteri  or  di  ameters  showed  a 
smaller  variation  with  maximal  diameters  observed  between  14  and  21  days  after  fracture.  These  gross 
measurements  can  be  used  to  develop  a  general  understanding  of  the  overall  rate  and  spatial  distribution  of 
tissue  formation  (peak  size  of  calluses)  and  cal  I  us  remodeling  (changein  size  and  progression  back  to  control 
values)  aid  provide  a  means  for  relating  quantifiable  histological  parameters  to  biomechaiicd  properties 
including  strength,  stiffness,  and  moments  of  inertia 

The  overall  progress  on  of  skeletal  morphogenesis  was  next  assessed  by  measurement  of  the  resorption 
of  cartilage  during  the  endochondral  phase  of  healing.  Total  osseous  tissue  (existing  cortical  bone,  new 
trabecula  bone,  new  osteoid  tissue,  and  the  void  space  that  includes  the  bone  marrow  cavity,  hematopoietic 
el  emails  and  empty ,  unstai  ned  space)  i  s  a  reason abl  e  measure  of  overal  I  pri  mary  and  secondary  bone  formati  on 
(Figure  3).  Measurement  of  void  areain  relationship  to  the  amounts  of  total  osseous  tissues  and  cartilagewas 
particularly  informative  as  this  measurement  provided  anotha  gross  assessment  of  the  rates  of  initial  cartilage 
resorption.  This  is  becausethe  void  spae,  which  is  composed  of  hematopoietic  elements  and  marrow  space, 
does  not  appear  withi  n  the  cal  I  us  unti  I  mi  neral  ized  carti  I  age  begi  ns  to  be  remodel  ed. 

The  most  striking  aspect  of  these  data  was  that  the  formation  of  the  carti  I  age  tissues  takes  place  in  an 
a^mmetric  manner  with  more  carti  I  age  formed  on  the  distal  than  the  proximo  side  of  thecal  us.  This 
asymmetry  was  also  reflected  in  the  overall  proximal  to  distal  area  and  diameter  measuements  of  thecdlusas 
seen  i  n  Fi  gure  2.  Compari  son  of  the  measurements  among  three  pri  mary  hi  stoi  ogi  cal  parameters  Cg,  TOT  aid 
Vd  at  14, 21,  and  35  days  demonstrated  that  spatial  asymmetry  in  the  intiaf  endochondral  distribution  effected 
both  the  pattern  of  tissje  remodeling,  and  the  pattern  of  osseous  tissue  formati  on  ova  the  time  course  of 
healing.  In  this  context,  the  measurement  of  void  aea  provided  a  retrospective  spatial  picture  of  whae 
remodel  i  ng  had  a  ready  occurred.  I  ntial  ly  the  void  aeas  wae  observed  a  the  sites  most  proxi  mal  anddista  to 
the  fracture.  Howeva ,  as  healing  preceded,  the  aeas  of  remodel  i  ng  progressed  i  nwad  and  the  voi  d  aeas  wae 
seen  flanki  ng  the  regions  adj  acent  to  the  fracture  at  the  carta  of  the  cal  I  us,  whae  the  aeas  of  cat  i  I  age  pasi  sted 
(Figure  3A  compaes  void  measurements  at  days  14  aid  35). 
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Us'ng  the  means  from  the  entrire  set  of  measurements,  the  percent  composition  of  cartilage  to  tots 
osseous  tisaje  to  void  area  were  determined  (Figure  3  panel  B).  While  it  has  generally  been  assumed  that 
cartilage  makes  14)  a  very  hi  gh  percentage  of  the  ti  ssue  content  of  fracture  callus  the  cart  i  I  venous  component 
comprised  only  ~10  to  15%  of  the  total  tissue  content  at  14  da/s  post  fracture.  However,  in  some  of  the  distal 
sections  doseto  the  center  cartilage  consituted  approximately  30%  of  the  cal  I  us  area  As  the  heating  progressed 
and  cartilage  was  resorbed  and  there  were  increases  in  the  percent  of  total  osseoustissue  and  void  area  This 
observation  regarding  thetotal  amount  of  cartilage  (10  to  15%)  isalso  consistent  with  that  seen  at  thepeefc 
period  of  cal  I  us  formation  in  mouse  femur  cal  I  uses  produced  in  the  same  experimental  maina  (data  not  shown). 

The  spat i a  pattern  of  tissue  remodeling  was  further  examined  by  assessing  osteoclast  density  (OC/Ar) 
throughout  the  callus  using  TRAP  staining  (Figure  3C).  Deta led  data  from  21  and  35  day  causes  a-e  presented. 
A  comparison  of  the  spatial  distribution  of  osteoclasts  with  the  spatial  distribution  of  cartilsgeard  void  aeasat 
21  da/s  provides  a  consistent  picture  of  the  retrospective  spatial  pattern  of  cartilage  tissue  resorption.  Theaeas 
containing  the  highest  dena'ty  of  osteodastswere  situated  toward  both  the  proximal  and  distal  edges  of  the 
cal  I  us  with  a  high  density  in  the  areas  of  carti  I  age  which  will  undergo  resorption.  As  these  areas  are  remodeled, 
th^  will  be  replaced  with  void  areas  and  trabecula  bone.  From  a  quantitative  paspedivethae  is  almost  a  5:1 
rati  0  i  n  the  total  densi ty  of  osteod  asts  duri  ng  the  adi  ve  pai  ods  of  carti  I  age  remodel  i  ng  at  day  21  than  duri  ng 
primary  and  secondary  bone  remodeling  that  predominates  at  day  35.  By  35  days,  the  distribution  of  osteodasts 
throughout  the  cal  I  us  is  relatively  uni  f am  except  at  the  most  proximal  edges. 

Three  Dimensional  Reconstruction  of  Callus  Morphogenesis 

Three  dimensional  reconst rudions  of  rat  fradure  calluses  wae  genaated  from  serial  histologies 
sections  taken  at  -100  Dm  increments  (all  60  sections),  spanning  -  6  mm,  of  a  14  day  callus  (Figure 4).  Four 
representative  si  ices,  two  from  the  proximal  and  two  from  the  distS  hal  ves  of  the  fradure  cal  I  us,  ae  presented 
in  the  bottom  panel  B.  The  three  dimensional  renderi  ngs  of  these  datavali  date  the  two  dimensions 
compositions  assessments  (see  below  Table  1)  and  more  dealy  show  that  the  catilage  components  of  the 
cS  I  us  wae  formed  in  an  asymmdric  manna,  with  the  majority  of  the  cartilage  famed  distS  to  the  fradure.  A 
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second  reconstruct  on  of  a  different  fourteen  day  cal  I  us  is  seen  i n  panel  C.  However,  thi s  renderi  ng  was 
stretched  so  that  a  1:5  to  1:0  aspect  ratio  is  presented  to  depict  the  asymmetry  of  the  tissue  in  a  more 
exaggerated  manner.  Two  separate  anatomical  perspectives,  viewing  thetissuefrom  either  the  medial  or  let  era 
surfaces,  again  show  the  asymmetry  in  cartilage  development  in  the  proximal  and  distal  orientations  but  edend 
this  the  asymmetry  into  the  medial /lateral,  posterior/anterior  dimensions  as  well.  The  cartilaginous  component 
of  the  cal  I  us  wraps  around  the  cortical  bone  toward  the  media  side  as  it  devd  ops  toward  the  proxima  end  of  the 
caius  Whereas  the  signdsthd  initiate  repar  ether  arisefrom  the  marrow  or  are  released  from  theinjured  bone 
marix  around  the  fracture,  thespaia  morphogenesis  of  the  caius  is  reguiaed  by  other  extrinsic  or  intrinsic 
factorstha  lead  to  asymmetry  of  tissue  development.  The  asymmetric  configurai  on  of  the  caius  continues  to 
be  observed  throughout  the  rest  of  the  period  of  the  bone  heaing  as  seen  in  the  three  dimensiona  renderings  of 
the  caius  a  21  cteys,  with  very  little  cartilage  now  remaning  on  the  proxima  side  of  the  fracture  (Figure  4D). 

An  assessment  of  how  thefracturecaiuseswerespaiaiy  remodeled  over  time  was  carried  out  by  three 
dimensiona  reconstructions  35  days  post  fracture  (Figure  5).  One  of  the  most  striking  feaures  of  these 
reconstructions  was  observed  in  thespaia  geomary  of  the  newly  formed  bone  within  the  35  day  caius.  By  35 
days  no  observable  cartilage  remaned  in  the  caius  (Figure  5A-B).  These  audies  show  tha  the  origina  centra 
space  of  the  caius,  which  was  occupied  by  cartilage,  has  now  been  replaced  with  trabecuiaed  bone  tissue  tha  is 
surrounded  by  an  outer  and  i  nner  she!  I  of  nav  bone  (Fi  gure  5C) .  The  outer  sha  I  formed  as  new  bone  grew  over 
the  cartilage  and  became  the  new  periostea  margin  adjacent  to  the  muscle.  The  inner  sha  I  formed  as  a  new 
surface  adjacent  to  the  origina  cortex  and  is  representaive  of  the  periostea  response  (Figure  5C).  The  origina 
cortices  tha  underlied  the  new  bone  surface  were  now  undergoing  extend  ve  resorption  and  becomethesiteof 
na/v  blood  vessa  ingrowth.  Eventudly,  both  the  origina  cortices  and  the  new  layer  of  bone  became  remodaed 
and  were  replaced  with  a  continuous  volume  of  trabecuiaed  bone  underlying  the  new  cortica  shai.  The 
asy mmary  of  the  bone  formai on  conti  nues  to  be  obvi  ous  duri  ng  this  phase  of  ca I  us  morphogenesi  s  as  seen  i  n 
the  cross  sectiond  views  There  is  a  very  smdl  amount  of  trabecuiaed  bone  bridging  the  inner  and  outer  sha  I  of 
the  cortices  i  n  the  proxima  view  but  an  extensive  network  of  trabecuiaed  bone  when  viewed  from  the  distd 
orientaion. 
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The  three  di  mensional  renderi  ngs  were  used  to  cal  cul  ate  the  percent  compositi  on  of  the  three  maj  or 
ti  ssue  types  within  the  cd  I  us.  These  were  based  on  comparing  the  determined  volumes  with  the  average  area 
fractions  calculated  from  the  two  dimensional  analyses  carri  ed  across  the  seri  al  sections.  The  results  for  these 
analysesaresummarized  in  Table  1  for  determinations  madefromthe21  d^  specimens  Ascanbeseen,  the2 
dimensional  analysis  accurately  (within  1-2%)  reflects  the  determined  percentages  of  3  dimensional  volumes 
val  i  dati  ng  both  the  i  ntemal  cona  stency  and  accuracy  of  both  sets  of  determi  nati  ons  made  from  the  hi  stol  ogi  cal 
measurements. 

Three  Dimensional  Reconstruction  of  the  Spatial  Pattern  of  Collagen  mRNA  Expression 

Using  in  situ  hybridization  the  spatial  pattern  of  collagen  types  I  and  II  gene  expression  were  three 
dimensionally  reconstructed  as  a  means  of  identifying  the  spatial  patterns  of  biological  activity  of  specific 
populati  ons  of  cells  in  the  cal  I  us  tissue.  For  these  studies  murine  femur  fractures  were  generated.  Both  light  and 
dark  field  images  show  the  fidelity  erf  the  in  situ  hybridizationsforthecollal  andcoi2a1  mRN As  (Figure  6). 
The  intense  and  very  specific  localize  on  of  the  two  probes  with  respect  to  either  bone  lining  cel  I  son  the 
nascent  trabecular  surfaces  or  at  the  interface  with  chondrocytes  adjacent  to  these  surfaces  are  shown. 

In  Figure  7,  murine  femur  fracture  cal  I  us  was  reconstructed  using  sections  taken  at  -500  Dm  increments 
spanning  ~6  mm  of  the  callus  Representative  proximal  and  distal  serial  sections  are  presented  (using  light 
microscopy/safranin  O  staining  for  orientation)  for  each  hybridization  probe  shown  in  Panel  A.  The  lower 
panels  depict  both  the  reconstruction  of  the  cortical  bone  in  region  to  cartilage  (Panel  B)  and  the  expression  of 
collAI  and  col2a1  inthesetissues(Panesl  D  and  E).  An  end  on  image  is  shown  to  provide  confirmation  of  the 
al  i  gnment  of  seri  al  reconstruct  ons  based  on  the  positi  on  of  the  medul  I  ary  pi  n  (Panel  C) .  M  ouse  femoral 
fracture  callus  showed  the  distal-lateral  asymmetric  distribution  of  cartilage  si  milar  to  that  observed  for  the  rat 
femur.  As  expected,  areasof  col2a1  expression  were  spatially  aligned  with  thecartilage.  The  expression  of 
col  2a1  was  confi  ned  to  the  central  areas  around  the  fracture  I  i  ne  but  not  the  areas  of  cal  I  us  more  peri  pherd  to 
the  fracture.  In  contrast,  collal  expression  was  uniformly  seen  on  either  side  of  the  fracture  and  was  distinctly 
excluded  from  the  regions  around  the  fracture  site.  Thecollal  expression  was  spatially  adjacent  to  ail  of  the 
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areas  in  which  carti  I  age  was  observed  aid  coil  al  was  expressed  i  n  a  mi  rror-l  i  ke  fashi  on  on  opposite  surfaces 
from  cartilage. 

DISCUSSION 

The  morphogenesis  of  skeletal  tissues  during  development  iscontrolled  by  spatio-temporal  relationships 
between  various  cellstypes  and  their  secreted  morphogens,  as  well  as  mechanical  cues  in  the  microenvironment. 
During  long  bone  development  involving  an  endochondral  process  these  spatial  arrangements  are  defined  by 
both  the  patterns  of  I  inear  growth  of  the  tissues  as  well  astheweight-bearingdemandsthesetissueswill  be 
required  to  meet.  While  it  is  generally  accepted  that  the  endochondral  bone  formation  during  fracture  healing 
recapitulates  many  of  the  molecular  events  that  control  embryological  and  post-natal  growth123,  the  well-defined 
spdiai  arrangements  thd  define  tissue  morphogenesis  during  bone  growth  are  not  as  easily  discernable  within  a 
fracture  cal  I  us.  The  question  that  arises  is  how  the  events  of  spatial  morphogenesi  s  intri  nsic  to  long-bone 
development  differ  from  those  thd  are  involved  during  healing  and  regeneration. 

Most  histological  assessments  of  fracture  cal  I  us  formation  have  used  longitudinal  sections  (sagittd , 
coronal)  whi  ch  have  I  ead  to  a  si  mpl  i  Stic  view  of  cal  I  us  morphogenesi  s  as  a  uniform  ri  ng  of  carti  I  age  that  i  rati  ally 
forms  around  the  fracture18.  However,  many  long  bones  are  ellipsoid  or  have  flattened  surfaces  aid  these 
configurations  ^fect  the  sp^ial  parameters  of  morphogenesis.  I  ndeed,  in  the  course  of  our  attempts  to  perform 
histomorphometric  analyses  of  fracture  healing,  we  observed  that  measurements  of  cartilage  aid  bone 
compositi  on,  as  well  as  measurements  of  cal  I  us  secti  onal  area  van  ed  tremendously  among  I  ongitudi  nai 
sections17  In  order  to  circumvent  these  problems,  a  serial,  transverse  sampling  method  was  used  in  the 
experiments  reported  in  this  study. 

The  moa  striking  finding  of  the  quantitative  analysis  of  tissue  composition  wasthe  degree  of  spatial 
variability  of  carti  I  age  formed  during  fracture  healing.  The  carti  I  age  tissue  spiraled  around  the  femoral  shaft  in  a 
proximal  to  distal  manner.  This  spatial  pattern  suggests  that  thereisa  much  greyer  complexity  among  the 
biological  and  biomechanical  factors  that  initiate  and  regulate  the  endochondral  processes  of  fracture  heeling 
compared  to  those  that  occur  in  the  later  phases  of  bone  remodeling.  This  inherent  complexity  is  reflected  in  the 
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much  greater  number  of  uniquely  expressed  mRNAs  in  fracture  callus  compared  with  unfractured  bone  as 
determined  by  large  scale  transcriptional  profiling.50  This  variability  m^/ also  be  the  cause  of  the  statistical 
variation  in  the  biomechanical  properties  of  endochondral  fracture  caHuset  early  time  points.2122-23  The  most 
likely  explanation  for  this  variation  is  related  to  the  anisotropy  produced  by  spatial  irregularities  in  the  formation 
of  cal  I  us  tissues.  Finally,  it  m  a/  be  speculated  that  the  complexity  of  theendochondra  process  makes  the 
transition  through  this  phase  of  bone  healing  an  extremely  sensitive  index  of  the  progression  of  normal  healing. 

I  ndeed,  delays  or  failures  to  progress  through  the  chondrogenic  stages  of  fracture  healing  ha/e  been  associated 
with  the  development  of  delayed  and  nonunions  in  both  clinical  studies18  and  investigations  in  animals17,24,25 
Two  of  the  most  interesting  findings  from  the  present  studies  are  the  reproducible  asymmetric  pattern 
i  n  w  hi  ch  carti  I  age  devel  ops  duri  ng  cal  I  us  f  ormati  on  and  the  spati  a!  mechani  smsbywhichthecallusis 
remodeled  back  to  its  pre-fracture  size.  In  the  case  of  the  spatial  pattern  of  carti  I  age  tissue  development,  our 
analyses  demonstrated  that  the  majority  of  the  carti  I  age  formed  distal  ly  and  laterally  to  the  fracture  with 
periosteal  bone  formation  being  the  greatest  on  proximal  lateral  surfaces.  This  asymmetry  was  conserved  across 
species  in  both  mice  and  rata  The  initial  asymmetry  in  tissue  formation  subsequently  affected  the  spatial 
pdtemsof  cartilage  resorption  that  were  observed  as  the  chondrocytes  unden/vent  hypertrophic  m^urati  on. 

Thus  regions  of  chondrocyte  maturation  and  subsequent  resorption  progressed  from  the  edges  inward  with 
regions  flanking  the  proximal  fracture  undergoing  resorption  first  followed  by  those  areas  distal  to  the  fracture. 
As  to  the  rol  e  of  i  ntri  nsi  c  versus  extri  nsi  c  factors  control  I  i  ng  the  uni  que  pattern  of  ti  ssue  f  ormati  on  in  the  cal  I  us 
it  is  important  to  note  that  femurs  of  several  mammalian  species  appear  to  grow  predominantly  from  the  distal 
physia26  27 28  Whi  le  we  do  not  report  on  data  from  the  ti  bi  a  we  have  observed  i  n  other  expert  merits  that  tibid 
fracture  cal  I  us  shows  more  carti  I  age  withi  n  the  proxi  mal  hal  ves  of  the  cal  I  us  and  this  bone  grows  preferentidly 
from  the  proximal  end.27  These  data  suggest  that  the  proximal  to  distal  distribution  of  cartilage  may  bean 
intrinsic  feature  of  the  way  skeletal  tissuesform  during  fracture  healing  and  msy  be  bone-specific  that  se 
related  to  additional  factors  that  are  developmental  controlled  such  as  blood  supply,  soft  tissue  or  muscular 
coverage 9  30  and  mechanical  demand.  Indeed,  oneofthemost  important  extrinsic  factors  is  the  biomechanical 
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environment,  which  has  been  shown  both  through  theoretical  cons' derations 31,32 and  empirical  experimentation 
to  be  crucial  to  cartilage  development. 33-37 

The  final  aspect  of  these  ddato  discuss  isthe  spatial  pattern  by  which  woven  bone  is  remodeled  to  its 
original  cortical  dimensions  Once  the  cartilage  component  of  the  cal  I  us  has  been  resorbed  there  are  two  new 
surfaces  of  bone:  an  i  nner  surface  that  has  grown  over  the  ori  gi  nal  cortex  and  an  outer  thi  nner  I  ayer  that  has 
encapsulated  the  cal  I  us  and  forms  the  new  interfacewith  the  periosteum.  A  trabecular  structure  bridges  these 
surfaces  i  n  the  space  that  was  occupi  ed  by  carti  I  age.  The  outer  she!  I ,  because  it  i  s  I  ocated  far  from  the  geometric 
center  of  the  bone,  will  be  responsible  for  the  majority  of  weight  bearing.38  This  outer  shell  is  connected  to  the 
original  cortex  via  trabecular- 1  ike  struts  which  appear  to  provide  sufficient  support  in  order  that  it  can  function 
to  Sabi  I  ize  the  fracture  This  represents  an  efficient  mechanism,  using  minimal  material,  to  rapidly  restore 
biomechanical  stiffness  and  strength,  while  allowing  for  remodeling  on  internal  surfaces.  A  similar  pattern  of 
morphological  expansion  has  been  observed  during  the  early  phase  of  rapidly  growing  animals  such  as  inbred 
mice39.  Unlike  long-bone  expansion  in  which  there  is  a  balance  between  periosteal  appositionai  growth  aid 
bone  resorption  at  the  endosteal  surface,  the  remodeling  of  fracture  callus  uses  different  and  unique  spatial 
mechanism  to  model  and  remodel  from  the  outer  surface  inwad  balancing  external  removal  with  the  addition  of 
bone  on  internal  surfaces.  Thisaspect  of  fracture  healing  is  perhapsthe  most  structurally  complicated  and  the 
ability  to  construct  a  simple  model  to  explain  how  this  takes  place  is  very  challenging.  As  noted  from  these 
results,  more  bone  forms  distal  to  the  fracture  site  during  femoral  healing  aid  proximal  to  the  fracture  site  in  the 
tibia  Further  investigations  will  be  required  in  order  to  better  understand  the  complex  interplay  of  each  of  these 
processes 
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FIGURE  LEGENDS 


Figure  l.  Schematic  presentation  of  methodological  standardization  for  histomorphometric  assessment  of 
fracture  cal  I  us.  In  these  black  and  white  pictures  areas  stained  red  with  the  Safranin  O-Fast  Green  appear  dark 
grey  to  black.  Pane!  A)  Diagram  of  the  sampling  approaches  used  to  quantify  histological  measurements  across 
a  fracture  cal  I  us.  The  fracture  cal  I  us  was  divided  into  segments,  which  were  denoted  as  proximd  or  distal  to  the 
center  of  the  callus.  Each  segment  was  100  pm  apart.  Sectionswere  sampled  within  segments  that  were 
progressively  spaced  at  larger  increments,  every  segment,  every  other  segment,  etc.  as  denoted  by  theXsacross 
the  whol  e  proxi  mal  to  di  sta!  ori  entati  on  of  the  cal  I  us  (M  atari  al  s  i  n  M  ethods,  Stati  sti  cal  M  ethods) .  Paiel  B) 
Representative  set  of  four  transverse  serial  histological  sections  taken  at  -2000  pm  increments  in  the  proximal 
(F)  to  di  stal  (D)  ori  entati  on  from  the  fracture  I  i  ne  of  a  14  day  post  fracture  cal  I  us  of  a  rat  femur. 
Photomicrographs  were  taken  at  25  x  magnification.  Panel  C)  Sampling  method  that  was  used  to  obtan 
mi crographs  for  analysis  for  measurements  of  osteod asts.  A  representati  ve  transverse  secti on  i s  presented. 
Micrographs  were  taken  consecutively  at  6-7  abitrary  positions  as  boxed  across  the  diameter  of  the  cal  I  us. 

T ransverse  i  ncrements  ae  then  taken  across  the  enti  re  length  of  the  cal  I  us  withi  n  segments  at  1 000  pm 
increments. 

Figure  2.  Compaison  of  mean  callus  aea  measurements  at  14, 21  and  35  days  post  fracture  in  rat  femur 
fracture  cal  I  uses.  Panel  A)  Graphical  analysis  of  the  mean  total  callus  aeas  from  1000pm  increments  from  the 
proximal  to diaal  sidesof  the  cal  I  us  as  a  f  unction  of  time  afta  fracture.  Panel  B)  Comparison  of  the  mean  cross 
sectional  measurements  of  callus  diameters  in  the  medial/lateral  or  posterior/anteri or  dimensions  as  a  function  of 
time  after  fracture.  Pane!  C)  Comparison  of  the  mean  cal  I  us  areas  as  a  function  of  time  after  fracture.  Vdues 
are  based  on  the  mean  taken  from  the  mi  ni  mal  i  ncrement  withi  n  the  95%  confi  dence  i  ntervsl  of  the  reference.  | 
p<  .05  between  35  days  and  14  days,  *  p<  .05  between  35  days  and  21  days. 

Figure  3.  Comparisons  of  the  histomorphometric  measurements  of  rat  femur  callustissue  compositions^  14, 

21,  and  35  days  post  fracture.  Panel  A)  Graphical  analysisof  the  primary  histological  measurements  assessed  at 
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500 rim  increments  from  the  proximal  to  distal  sides  of  the  callus.  Panel  B)  Graphical  analysis  of  Cg,  TOT  aid 
Vd  as  a  function  of  time.  Comparison  of  the  mean  cal  I  us  areas  as  function  of  time  after  fracture.  vauesae 
based  on  the  mean  taken  from  the  mi  ni  mal  i  ncrement  withi  n  the  95%  confidence  i  nterva  of  the  reference.  Error 
bar  =  SD  of  the  means  from  three  cal  I  uses  at  each  ti  me  point  t  P<  05  between  35  days  aid  14  days.*  p<  .05 
between  35  days  and  21  d^s.  Panel  C)  Quantification  of  osteoclasts  within  fracture  cal  I  us  at  21  d^spost 
fracture.  Left  panel  shows  the  mean  total  numbers  of  osteocl  asts  per  area  at  1 000  pm  i  ncrements  across  the 
proximal  to  distal  length  of  the  cal  I  uses  from  the  center.  Values  from  ather  21  or  35  days  post  fracture  ae 
shown.  Right  panel  presents  overall  mean  osteoclast  number  pa  unit  aea  when  avaaged  ova  the  entire  length 
of  the  cal  I  us.  I  ncrements  ae  avaaged  val  ues  f  rom  the  mean  additi  ve  val  ues  from  the  proxi  mal  and  di  stal 
increments  at  the  centa  1000, 2000,  and  3000  pm.  Error  ba  =  SD. 

Figure  4.  Reconstruction  of  a  rat  femoral  fracture  callus,  14  days  post  fracture.  Three  dimension^ 
reconstructions  wae  made  from  saial  sections  taken  every  100  microns  Panel  A)  Representative  3D 
reconstruction  of  thecartilagetisajeswithinthecallusrelativetotheoriginal  cortical  bona  Thecortica  boneis 
denoted  by  grey  to  white  staned  aeas  relative  to  mature  cartilage  denoted  dak  grey  to  black.  The  positions  of 
the  levels  of  individual  slices  depicted  in  panel  A  are  corraaed  with  those  in  Pare!  B)  The  aiaomica 
orientaionsof  four  ^ices  through  the  cal  I  us  tissue.  The  positions  of  the  a  ices  as  denoted  in  panel  A  ae  each 
shown  in  panel  B.  Panels  C)  Demonaration  of  the  geometric  asymmetry  of  the  cal  I  us  formation.  Media  aid 
laieral  presentations  of  a  second  reconstruction  of  a  14  day  rat  cailusae  shown  in  a  1:1.5  aspect  ratio.  Arrows 
denote  orientation  of  rotation  and  volume  distribution  of  cartilage.  Panel  D)  Reconstruction  of  the  cortica  bone 
(light  grey  to  white)  raaivetomaurecartilage(dak  grey  to  black)  withinthesoft  caiusof  a21  day  caius 
shown  for  femur  fracture  of  the  ra. 

Figures.  Reconstructions  of  ra  fracture  ca  I  uses  a  35  days  Panel  A)  represents  cortica  bone  as  described  in 
figure  6.  No  cartilage  was  observed  in  this  specimen.  Panel  B)  SI  ices  from  selected  proxi  ma  aid  distal 
positions  are  as  denoted  in  figure  6.  Panel  C)  Reconstructions  of  ra  femur  caiuses  a  35  days.  The  left  paid 
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depicts  an  end  on  view  from  the  distal  orientation  of  the  entire  cal  I  us  volume.  The  original  corticd  bone  is 
denoted  with  (*).  Middle  and  right  panels  present  proximal  and  distal  views  of  the  external  shell  of  new  bone 
and  the  trabecular  structure  that  has  formed  during  the  resorption  of  the  regions  of  original  endochondral  bone. 
Boxed  areas  are  regions  in  which  secondary  remodeling  has  completely  replaced  the  original  corticd  bone  with 
a  continuum  of  trabecular  bone. 

Figure  6.  Three  dimensional  reconstruction  of  mouse  femur  fractures  and  collagen  types  I  and  II  gene 
expresaon.  Panel  A)  Represents  ve  light  and  dark  field  images  of  micrographs  taken  of  in  situ  hybridizations. 
Top  micrograph  depicts4Gx  magnification  light  field  of  a  Saffrin  O  /  Fast  Green  staned  section  from  a  7  day 
muri  ne  femur  fracture  cal  I  us.  M  iddl  e  and  bottom  panel  s  are  dark  f  i  elds  from  seri  al  secti  ons  hybri  dized 
respectively  with  col  lal  andcol2a1  DNA  probes.  Boxed  area  corresponds  to  that  depicted  in  the  lOOx 
magnifies  ons  seen  in  panel  son  the  Right.  AlternSng  images  from  top  to  bottom  show  light  field  and  dale 
fields  from  consecutive  serial  si  ides  hybri  dized  respectively  with  a  probe  for  col  lal  and  col2a1.  Identical  aeas 
of  silver  grain  development  seen  in  the  light  field  and  matched  dak  fields  ae  indicated  with  arows. 

Figure  7.  Spatial  reconstruction  of  tis&ie  structure  and  collagen  gene  expression  of  7  day  murine  femur  fracture 
calluses  Panel  A)  Micrographic  images  (4Qx  magnification)  of  light  (saffrin  O /  fast  green)  aid  dak  field  (in 
situ  hybridization)  images  of  saial  sections  of  calluses  made  at  500p.M  increments  ova  a  6  mm  distaice. 
Representative  saial  histological  wae  hybridized  with  probes  for  col  lal  aid  col2a1.  The  said  sections  wae 
then  reconstructed  Panel  B)  Reconstruction  of  the  spatial  relationship  of  thecartilage  tissue  within  the  cdlus  in 
relation  to  the  cortical  bone.  Panel  C)  Flat  surface  reconstruct! on  showing  end-on  view  of  thetotd  volume. 
Panel  D)  col  lal  and  col2a1  mRNA  expression  as  determined  by  in  situ  hybridization  in  relation  to  the  corticd 
bone.  Panel  E)  Thereldionshipof  collal  mRNA  expression  relative  to  cartilage.  Corticd  bone  is  cola  coded 
yellow,  cartilage  red,  aeas  of  col  lal  expression  blue  and  col2a1  black.  The  line  of  the  fracture  is  depicted  by", 
*  denertes  the  hole  from  the  intramedullary  pin. 
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Tabic  1.  Comparison  of  Area  Composition  Versus  Volume  Compositions  in  21-Day  Rat  Fracture  Calluses 


Sample 

Cartilage 

Callus 

Cortical  Bone 

Space 

1 

16.819 

169.084 

35.154 

10.712 

2 

20.249 

183.589 

27.109 

18.613 

3 

10.832 

158.645 

20.735 

10.868 

Mean 

15.967 

170.439 

27.666 

13.398 

Standard  Deviation 

4.766 

12.527 

7.226 

4.517 

Percentages 

9 

16  (83*) 

8 

Percentage# 

Histomorphometry 

9.0 

82 

9.0 

+  Rounded  to  the  nearest  whole  number 

*  Total  percentage  Cortical  bone  pi  us  all  other  Non  Cartilage  Osseous  Tissues 

#  H  i  stomorphometry  measurements  taken  from  mean  values  of  serial  sections  taken  at  every  500nm  over  the 
same  proximal  to  dteal  dimensions  as  that  which  was  reconstructed.  Hi  stomorphometry  measurements  end 
tissue  vol  umes  were  obtained  by  separate  methods  of  digital  irrcge  entry. 
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